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1. Introduction 

Hadron physics, ranging from particle to nuclear physics, aims to a precise and 
consistent description of hadronic structure and interactions. This is a formidable task 
that can only be ideally accomplished within the framework of the Standard Model 
of the strong and electroweak interactions. The strong interaction part, Quantum 
Chromodynamics (QCD), should give account of the main bulk of hadronic data but 
in its current state of development this objective seems far from being attainable. 
This has motivated the formulation of alternative descriptions, based on QCD and/or 
phenomenology, which restrict their study to a particular set of hadronic systems 
and are able to reproduce the data and to make useful testable predictions. In this 
sense non-relativistic constituent quark models incorporating chiral potentials provide 
undoubtedly the most successful, consistent and universal microscopic description of 
the baryon spectrum and the baryon-baryon interaction altogether. Therefore they 
are ideal frameworks to analyze few-baryon systems. The purpose of this review article 
is to report the progress made in the last years in this analysis. But first to put these 
models in perspective we shall start by recalling the main aspects of the different 
approaches in hadron physics. Then we shall present the historical development of 
chiral constituent quark models emphasizing their connection to phenomenology and 
their plausible relation to the basic theory. 

QCD [1] is nowadays accepted as the basic theory to describe the hadrons and 
their strong interactions. QCD is a renormalizable quantum field theory of quarks 
and gluons based on a local gauge principle on the group SU(3)coiour- Apart from 
this local symmetry the QCD Lagrangian has also, for massless quarks, a global 
SU(n) x SU(n) (n: number of quark flavours) chiral symmetry that appears to 
be spontaneously broken in nature. Other relevant properties of QCD are related 
to the running of the quark-gluon coupling constant: for high momenta (i.e., high 
momentum transfers as compared to the QCD scale, Aqcd) the coupling tends to 
vanish and the quarks are essentially free, a property named asymptotic freedom; for 
low momenta the coupling becomes strong, a property known as infrared slavery. 
As a consequence of this behaviour a solution of the theory is only attainable, 
perturbatively, at high momenta. In the low-momentum regime one has to resort 
to non-perturbative calculation methods (sum rules, instanton based calculations,...) 
of limited applicability or to reformulations of the theory from which to approach the 
exact solution. Thus lattice QCD has been developed in the hope of getting from it the 
exact solution by a limiting procedure from the discrete lattice space to the continuum. 
However progress in making precise, detailed predictions of the physical states of the 
theory is slow (due in part to the enormous computer capacity required) [2] . 

Alternatively effective field theory approaches have been proposed: from the 
action and measure of QCD one can integrate out the irrelevant degrees of freedom 
in the momentum region under consideration to obtain a more tractable field theory 
with the same S matrix. This procedure, used with notable precision at high momenta 
where the smallness of the coupling constant makes feasible to calculate the effective 
Lagrangian, cannot be pursued with the same precision at low momenta. Such 
a difficulty can be overcome by the application of a celebrated, though unproven, 
theorem by Weinberg [3] which states that one should write the most general 
Lagrangian, constructed from the accessible degrees of freedom in the momentum 
region under consideration, which satisfies the relevant symmetries of the theory. This 
is the scheme used in Chiral Perturbation Theory [4] which has been successful in the 
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study of low momentum hadron (in particular meson) physics. 

Other approaches take limits of QCD and generate from them systematic 
expansions to get corrections. Among them we shall briefly comment on the Heavy 
Quark Effective Theory (HQET) [5], the Non-Relativistic QCD (NRQCD) [6] and the 
1/N C approach [7,8]. 

For physical systems with a heavy quark, interacting with light quarks and gluons 
carrying a momentum of order Aqcd, an appropriate limit of QCD involves taking 
the heavy quark mass going to infinity, toq — > oo (HQET). In this limit the strong 
interactions of the heavy quark become independent of its mass and spin. These heavy 
flavour and spin symmetries, not present in QCD, lead to model independent relations 
that allow for instance the description of exclusive decays in terms of a few parameters. 
Corrections are obtained from an expansion in terms of 1/toq, the inverse of the heavy 
quark mass. 

When more than one heavy quark is present in the system (for example in 
quarkonia) the heavy quark kinetic energy, treated as a small correction in HQET, 
cannot be treated as a perturbation. For such systems the appropriate limit of QCD to 
examine is the c — > oo limit (or v/c — > in quarkonia, being v the relative heavy quark 
velocity and c the speed of light). The resulting field theory is called Non-Relativistic 
QCD. An expansion in terms of 1/c is done. NRQCD has improved the understanding 
of the decays and production of quarkonia. 

The 1/N C approach is the limit of QCD when the number of colours N c tends 
to infinity. An expansion in terms of l/N c , the inverse of the number of colours, is 
performed. Remarkably, large N c QCD reduces smoothly to an effective field theory 
of non-interacting mesons to leading order. Baryons emerge as solitons in this weakly 
coupled phase of mesons. This is the base of topological models where baryons are 
constructed from non-linear meson fields in chiral Lagrangians (Skyrme model [9], 
chiral soliton model [10],...) 

At a much more phenomenological level model building has been extremely useful 
to provide a classification of hadronic data and to make predictions. Though the 
connection with QCD is not clearly established and there is not a sound systematics 
to obtain corrections, models provide simple physical pictures which connect the 
phenomenological regularities observed in the hadron data with the underlying 
structure. 

Historically the first model of hadron structure, a non-relativistic quark model, 
appeared in the nineteen sixties right after Gell-Mann and Zweig [11,12] introduced 
the quarks as the members of the triplet representation of SU (3) Flavour from which the 
hadron representations could be constructed. The lack of observation of free quarks in 
nature seemed to indicate their strong binding inside hadrons. Dalitz [13] performed a 
quark model analysis of hadronic data based on flavour. This analysis seemed to point 
out some possible inconsistencies with the Spin-Statistic theorem. The introduction 
of colour [14, 15], a new degree of freedom for the quarks, saved the consistency and 
led to the formulation of QCD. The other way around QCD provided from its infrared 
slavery property some justification for coloured quark confinement inside colour singlet 
hadrons. Moreover an interquark colour potential was derived by de Riijula, Georgi 
and Glashow [16] from the one-gluon exchange (OGE) diagram in QCD, and applied 
in an effective manner to reproduce energy splittings in hadron spectroscopy. 

In the 70's the potentialities of quark model calculations in hadronic physics were 
established. The non-relativistic quark model was formulated under the assumption 
that the hadrons are colour singlet non relativistic bound states of constituent quarks 
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with phenomenological effective masses and interactions. 

Non-relativistic quark models incorporating effective OGE and confinement 
potentials provided in the heavy meson sector a very good description of the 
charmonium and bottomonium spectra [17,18]. The reasonable results obtained when 
applied to the light baryon sector [19] were somewhat surprising since according to 
the size of the baryons (~ 0.8 fm) and the mass of the light constituent quarks 
employed (m u — rrid — 300 MeV) the quark velocities were close to c and the non- 
relativistic treatment could hardly be justified. In the spirit of quark potential model 
calculations this pointed out that the effectiveness of the parameters of the potential 
could be taking, at least partially into account in an implicit form, some relativistic 
corrections. Years later the low-lying states of the meson [20] and baryon spectra [21] 
were reasonably reproduced with a unique set of parameter values. 

On the same line relativistic quark models started to be developed. The bag 
model [22] considered asymptotically free quarks confined into a space-time tube 
through field boundary conditions. Chiral symmetry was incorporated through current 
conservation at the boundary associated to the presence of a meson cloud [23] . Hadron 
spectra and properties were calculated. However, despite the very intuitive image 
the bag model provided of a nucleon as a quark core surrounded by a meson cloud, 
calculations had to face the technical difficulty of separating the center of mass motion, 
an endemic problem in all relativistic quark models. 

Mid-way between non-relativistic and relativistic treatments semirelativistic 
quark models, where the relativistic expression for the kinetic energy was used in 
the Schrodingcr equation, were explored [24]. 

At the same time there were attempts to study the hadron-hadron interaction 
within the non-relativistic quark model framework [25]. It was soon realized that the 
colour structure of the OGE interquark potential could provide an explanation for the 
nucleon-nucleon (NN) short-range repulsion [26]. This encouraged several groups to 
undertake the ambitious project of describing the NN interaction from the quark-quark 
(confincmcnt+OGE) potential [27-33] by using the resonating group method or Born- 
Oppenheimer techniques. According to the expectations a quantitative explanation 
for the short-range NN repulsion was found. However, for the medium and long-range 
parts of the NN interaction, though (qq) and (qq) 2 excitations generated by off-shell 
terms of the Fermi-Breit piece of the OGE were considered so that the NN potential 
became attractive in the 0.8 — 1.5 fm range, the attraction was too weak to bind the 
deutcron or to fit the extreme low-energy S'-wave scattering [34] . 

To remedy this hybrid quark models [35, 36] containing both, interquark 
(confinement+OGE) and interbaryon long-range one-pion and medium range one- 
sigma, or two-pion exchange potentials, were introduced. Although these models were 
efficient to reproduce scattering and bound state data it was compelling, for the sake 
of consistency, to find a justification for these interbaryon potentials at the quark 
level. By considering that they could be due to (chiral) meson clouds surrounding 
the nucleons quark core, and pursuing the philosophy initiated in [16] of trying to 
incorporate into the quark potential model the dynamics and symmetries of QCD, an 
implementation of chiral symmetry at the quark potential level was needed. 

To accomplish this task the progress in the understanding of the connection 
between quark potential models and the basic theory was of great help. Manohar 
and Georgi [37] argued that the scale associated to confinement in a hadron (Aqon ~ 
100—300 MeV) being smaller than the one associated to chiral symmetry breaking 
(A <~ 1 GeV) would drive to a picture where quarks, gluons and pions coexisted in a 
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region of momentum. Using an effective field theory approach they obtained a model 
where constituent quarks (with a mass generated through chiral symmetry breaking) 
and gluons interacted via conventional colour couplings while quarks and pions did 
via a non-linear sigma model. By the same time Diakonov, Petrov and Yu got similar 
results concerning chiral symmetry breaking from a picture of the QCD vacuum as a 
dilute gas of instantons [38,39]. 

The effect of incorporating to the simple (confmement+OGE) model a one-pion 
exchange (OPE) potential at the quark level started to be analyzed. It was realized 
that the OPE potential gave rise to a NN short-range repulsion [40] to be added to 
the OGE one. As a consequence the quark-gluon coupling constant needed to fit the 
data got nicely reduced from its former effective value (>1) to a value (<1) much 
more according to the perturbative derivation of the OGE from QCD. The same 
conclusion came out from the OPE contribution to the A— N mass difference [41]. On 
the other hand the OPE potential generated, at large internucleonic distances, a NN 
conventional pion exchange interaction as needed. Nonetheless at the medium range 
it did not provide enough NN attraction to fit the data. 

The introduction of sigma as well as pion exchanges between the quarks, in 
the form dictated by a Nambu-Goldstone realization of chiral symmetry within the 
linear sigma model, allowed to overcome this problem. Thus, the first non-relativistic 
constituent quark model of the NN interaction incorporating a chiral SU(2) x SU(2) 
potential came out [42,43]. NN phase shifts (with the exception of P- waves) and 
deuteron properties were satisfactorily described. Furthermore, reasonable, though 
not precise, hadron spectra were predicted with the very same model (i.e. with the 
same set of parameters fitted from the NN interaction) [42,44]. This model has been 
usually referred in the literature as chiral constituent quark model (CCQM). We shall 
maintain this logo hereforth. 

Later on a SU(3) x SU(3) chiral quark model where constituent quarks interact 
only through pseudoscalar Goldstone bosons (GBE) was developed to describe very 
successfully the baryon spectra [45]. However its application to the NN interaction 
revealed a too strong tensor force, generated from the 7r-exchange, and the absence 
of the necessary medium-range attraction [46] . When properly implemented with the 
one-gluon exchange and scalar 5/7(3) x SU(3) Goldstone boson interactions, these 
models were successfully applied to the NN and nucleon-hyperon interactions [47,48]. 

The significant role played by the semirelativistic kinematics in the GBE model to 
fit the spectroscopy encouraged the implementation of a relativistic treatment in the 
CCQM. When done [49] a pretty nice description of the non-strange baryon spectra 
was obtained as well. 

The success in the description of the non-strange baryon spectroscopy and the NN 
interaction in a consistent manner makes the CCQM to be a powerful tool to treat, 
in a parameter-free way and on the same footing, other baryon-baryon interactions. 
This consistent treatment is mandatory in the study of few-baryon systems due to 
the intertwined role of nuclcons and resonances in them. In this article we review the 
application of non-relativistic constituent quark models to obtain the energy spectra of 
hadrons and the baryon-baryon effective potentials (we will not consider baryon form- 
factors or amplitudes for weak, electromagnetic or strong decays since they are very 
sensitive to relativity [50]). Though most of the results will refer to the CCQM, results 
from other models will be included for completeness. The order of the presentation is 
the following. In section|5]the derivation of the chiral part of the quark-quark potential 
from an effective Lagrangian, incorporating spontaneous chiral symmetry breaking, is 
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detailed. A one-gluon exchange potential, giving account of the residual (perturbative) 
colour interaction, completes the quark-quark potential whose parameters are fitted 
from the NN interaction and the non-strange baryon spectrum. In order to be able 
to generate a baryon-baryon potential from the quark-quark one and to deepen the 
understanding of the role played by quark antisymmetry, a variational two-baryon 
wave function is introduced in section A detailed explanation of the NN short- 
range repulsion based on the interplay between quark antisymmetry and dynamics is 
presented. The baryon-baryon potential for NN, NA, AA and NN*(1440) channels 
are constructed in section 0] and section |SJ Comparison to experimental data, when 
available, and to other model results, are shown. A coupled channel calculation for 
the NN system above the pion threshold is reported. In section HJ a thorough study of 
the non-strange baryon spectrum is carried out. Some related comments on the first 
experimental candidate for an exotic baryon, the + resonance, are also included. 
Consistency of the baryon spectrum wave functions with the variational baryon-baryon 
wave functions is shown. Section is devoted to the search for unstable two- and 
three-baryon resonances. A more complete calculation for the triton binding energy 
is presented. Finally in section |H1 we summarize the main results and conclusions. 
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2. The chiral constituent quark model. 

2.1. The quark- quark potential 

In non-relativistic quark models quark colour confinement inside colour singlet hadrons 
is taken for granted. Though confinement has not been rigorously derived from 
QCD, lattice calculations show in the so-called quenched approximation (only valence 
quarks) an interquark potential linearly increasing with the interquark distance [51]. 
This potential can be physically interpreted in a picture in which the quarks are 
linked with a onc-dimensional colour flux tube or string and hence the potential is 
proportional to the distance between the quarks, 

VcoN(fy) = ~ a c K ■ %j r ijt (2.1) 

where a c is the confinement strength, the A's are the SU(3) colour matrices, and the 
colour structure prevents from having confining interaction between colour singlets 
[52]. Hadron sizes correspond to a scale of confinement Acon ~ 100—300 MeV. On 
the other hand low-lying hadron masses are much bigger than light quark (up and 
down) current masses in QCD. So when dealing with hadrons one can reasonably 
assume as a good approximation the light quarks to be massless. 

In the limit of zero light-quark masses the QCD Lagrangian is invariant under the 
chiral transformation SU (2) l <8> SU {2)r. This symmetry would imply the existence of 
chiral partners, that is, for each low-lying hadron there would exist another one with 
equal mass and opposite parity || . This is not observed in nature % what points out to a 
spontaneous chiral symmetry breaking in QCD. As a consequence, the current quarks 
get dressed becoming constituent quarks and Goldstone bosons are generated [55]. 
Would the whole process be exact, one would end up with massless Goldstone bosons 
exchanged between the constituent quarks. In the real world chiral symmetry is only 
an approximate symmetry so one ends up with low-mass bosons (with masses related 
to the quark masses) exchanged between the constituents. 

The picture of the QCD vacuum as a dilute medium of instantons [38,39] explains 
nicely the spontaneous chiral symmetry breaking, which is the most important non- 
perturbative phenomenon for hadron structure at low momenta. Quarks interact with 
fermionic zero modes of the individual instantons in the medium and therefore the 
propagator of a light quark gets modified and quarks acquire a momentum dependent 
mass which drops to zero for momenta higher than the inverse of the average instanton 
size ~p. The momentum dependent quark mass acts as a natural cutoff of the theory. In 
the domain of momenta q < 1/p, a simple chiral invariant Lagrangian can be derived 
as [38] 

L = iHi^d^ - MU j5 )ip (2.2) 

where C/ 75 = exp(i7T • T75 / f v ) . n denotes a Goldstone pseudoscalar field, f v is the pion 
decay constant, and M is the constituent quark mass. The momentum dependence of 

|| There are two ways in which a symmetry of a Lagrangian manifests itself in nature [53]. The first 
one is the standard Wigner-Weyl realization when the generators of the symmetry group annihilate 
the vacuum. In this case nature exhibits the symmetry in the form of degenerate multiplets. The 
second one is the Goldstone realization and it corresponds to the case of a vacuum that is not 
annihilated by all the generators of the group. The symmetry of the Lagrangian is not evident in 
nature, one says it is hidden, and this is referred to as the spontaneous breaking of the symmetry. It 
is essentially the case for QCD [53,54]. 

^ The splitting between the vector p and the axial ai mesons is about 500 MeV (2/3 of the p mass) 
and the splitting between the nucleon and its chiral partner is even larger (940 — 1535) MeV. 
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the constituent quark mass can be obtained from the theory. It has been effectively 
parametrized as M(q 2 ) = m q F(q 2 ) [41,43] with 



F(q 2 ) 



A 2 



A 2 



(2.3) 



where A is the chiral symmetry breaking scale, A ~ 1/p and m q = M(q 2 = 0) sw 350 
MeV. We shall call hereforth m q the constituent quark mass. C/ 75 can be expanded 
in terms of boson fields as, 

1 



IP* = 1 + — r 



(2.4) 
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The first term generates the constituent quark mass and the second gives rise to 
a pseudoscalar (PS) one-boson exchange interaction between quarks. The main 
contribution of the third term comes from the two-pion exchange which is usually 
simulated by means of a scalar (S) exchange. From the non-relativistic approximation 
of the Lagrangian one can generate in the static approximation the quark-meson 
exchange potentials: 
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(2.6) 



where the i and j indices are associated with i and j quarks respectively, stands 
for the interquark distance, a c h is the chiral coupling constant, the ct's (t's) are the 
spin (isospin) quark Pauli matrices, tops and tos are the masses of the pseudoscalar 
and scalar Goldstone bosons, respectively. Sij = 3 (o^ • ■ r^) — <?j • dj is the 

quark tensor operator and Y(x) and H (x) are the standard Yukawa functions defined 
by Y(x) = e~ x /x and H{x) = (1 + 3/x + 3/x 2 )Y(x). 

For q 2 ^> A 2 one expects QCD perturbative effects playing a role. They mimic the 
gluon fluctuations around the instanton vacuum and are taken into account through 
the OGE potential [16]. From the non-relativistic reduction of the one-gluon-exchange 
diagram in QCD for point-like quarks one gets 



■ a s At • A, 
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2^ ^ " 
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3 — fjj/ro 
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4to 2 r 3 - y 



(2.7) 



where a s is an effective strong coupling constant. Let us realize that the contact 
term involving a Dirac 5(f) that comes out in the deduction of the potential has been 
regularized in the form 

1 e - r / r a 

*(f) -> — - (2-8) 



47rr 2 



where ro is a regularization parameter giving rise to the second term of (|2.7|) . This 
avoids to get an unbound baryon spectrum when solving the Schrodinger equation [56] . 
Thus the quark-quark interaction has the form, 



V qq(fij) = V C ON(n 3 ) + VoGE(^) + Vps(nj) + Vs(fij) . 



(2.9) 
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Such a model has an immediate physical interpretation. In the intermediate region, 
between the scale at which the chiral flavour symmetry is spontaneously broken, A ~ 
0.8 GeV, and the confinement scale, A CO n ~ 0.2 GeV, QCD is formulated in terms 
of an effective theory of constituent chiral quarks interacting through the Goldstone 
modes associated with the spontaneous breaking of chiral symmetry. For q 2 3> A 2 
gluon exchange is also important. 

It is worthwhile to note that vector meson-exchange potentials (p, lo) are not 
considered. The problem of unifying the quark exchange and meson exchange in 
the nuclear force has been a matter of discussion [57]. It has been shown that the 
pseudoscalar (it) and scalar (a) meson-exchange terms can be simply added to the 
quark-exchange terms without risk of double counting. However, the vector-meson 
exchanges (p, ui), which play an important role in meson-exchange models at the 
baryon level need some care. In baryonic one-boson exchange models, the w-meson 
provides the short-range repulsion of the NN interaction. In chiral constituent quark 
models this task is taken over by the antisymmetrization effects on the pseudoscalar 
exchange combined with OGE. Besides, the p-meson reduces the strength of the 
tensor pseudoscalar interaction, the same effect that is obtained from the quark- 
exchange terms of the pseudoscalar potential as has been checked in charge-exchange 
reactions [58]. 



2. 2. Model parameters 

In the spirit of quark model calculations the parameters in the potentials have an 
effective character. A rough estimate of the values of the parameters can be made 
based on general arguments. It is well established that the NN interaction at long- 
range is governed by the one-pion exchange. Therefore, to reproduce accurately this 
piece of the NN interaction, one is forced to identify the mass of the pseudoscalar field 
with the physical pion mass. The mass of the scalar field, the one-sigma exchange 
(OSE), is obtained by the PCAC relation [59] 

m 2 PS = ml (2.10) 
TOg = TOp S + 4m 2 , . 

As the pseudoscalar field is identified with the pion, the a c h coupling constant should 
reproduce the long-range OPE interaction. If two nucleons are separated enough, 
the central part of Vp S , the pseudoscalar interaction between quarks, generates an 
interaction between nucleons of the form, 

i A2 „-m„_R / r \ 2 

V PS (R) = -a ch G A2 - i^-j (a N -a N )(r N -r N ) (2.11) 

where R is the interbaryon distance and G depends on the nucleon wave function. 
Comparing with the standard one-pion-exchange internucleon potential, 

I f2 A2 p -m„R 

VoMR) = 3 ^ A2 _ ml R {Sn ■ a N )(r N ■ r N ) (2.12) 

where the same form factor has been used at the quark and baryon levels, and using a 
harmonic oscillator wave function for the nucleon in terms of quarks (see next section) , 
one finally obtains [60], 
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This gives the chiral coupling constant a c h in terms of the 7rNN coupling constant, 
taken to be /^ NN /4"7r = 0.078 [61], and b. Usual values in the literature for b range 
between 0.4 and 0.6 fm. The most stringent determination, b — 0.518 fm, was done 
in reference [62] by means of a simultaneous study of S-w&ve NN phase shifts and 
deuteron properties (see section |SJ|. As it will be discussed, this value turns out to be 
consistent with the solution of the baryon spectra as a three-body problem in terms 
of the interaction lj2.9jl [44]. The value of a c h is then a ch =0.0269. 

The tensor force is mainly due to the pion interaction. Therefore, the value of 
A can be fixed examining a process dominated by the tensor term. Such a reaction 
could be the pp — > nA ++ because, at high momenta, more than 90 % of the interaction 
corresponds to the tensor part. The calculation of reference [58] suggests for A a value 
close to 4.2 fm _1 . 

The value of a s is estimated from the A— N mass difference. In the chiral 
constituent quark model there are contributions not only from spin-spin term of the 
OGE but also from the pseudoscalar interaction, the latter contributing approximately 
half of the total mass difference. The rest is attributed to the OGE, and the value 
of a B is adjusted to reproduce the experimental A— N mass difference. This gives 
a s ~ 0.4 — 0.5 for tq — 0.25 fm, a standard value of ro within the stability region for 
the A— N mass difference (see section Ifl. 2 .4JI . 

The constituent quark mass m q is an effective parameter whose value is 
conventionally chosen in the range of 300 MeV, using in general smaller values for 
the study of the NN interaction and bigger ones for one-baryon properties. From 
the proton and neutron magnetic moments in the impulse approximation one gets 
uiq = 350 MeV. From NN scattering the quark mass value should be close to one 
third of the nucleon mass m q ~ 1/3tojv ~ 313 MeV. In this way a correct relation 
between the momentum and kinetic energy is guaranteed [27]. 

Finally, regarding a c we should first notice that the contribution of confinement 
to the force between two baryons is very small (zero for quadratic confining, see 
section I3.2JI . Hence its value is only constrained by the requirement of having a 
confining (a c > 0) and not a deconfining (a c < 0) interaction. This can be guaranteed 
through the nucleon mass stability condition dM^(b)/db\ b=Q 51g = 0. Concerning its 
specific value we can resort to the baryon spectrum which is strongly dependent on 
it (see section EJl. We quote in tabled the standard value derived from the baryon 
spectrum analysis. Some caution is necessary when comparing the strength of the 
confining potential to other values given in the literature. First, one should be aware 
of the specific form used for the confining interaction, if the colour Gell-Mann matrices 
are used or not (a factor —8/3 is in the way). Second, when scalar potentials between 
quarks are used, as it is the case of the chiral constituent quark model, smaller values 
of a c than in pure OGE models are obtained. 

The standard values of the parameters used in the chiral constituent quark model 
is resumed in tabled 



Table 1. Quark model parameters. 
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a s a c 
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3. The non-strange two-baryon system 

3.1. The two-baryon wave function: quark Pauli blocking 

The calculation of the interaction between two baryons requires the knowledge of the 
two-baryon wave function and therefore that of a single baryon in terms of quarks. 
Single baryon wave functions have been calculated using different methods available 
in the literature to solve the three-body problem (see section [fjj). Although the 
resulting wave functions have an involved structure, it has been shown that for the 
baryon-baryon interaction they can be very well approximated by harmonic oscillator 
eigenfunctions [44]. For the non-strange baryons we are going to consider: N, A and 
N*(1440), the total wave function of a single baryon can be explicitly written as the 
product of three wave function components: spatial, spin-isospin, and colour space 
respectively, + 

\$B) = (l>(f>i,r2,r3;R)®xsT<8>&, (3.1) 

where f\ is the position of quark i and R denotes the center of mass coordinate of the 
baryon. The symmetrization postulate requires the wave function to be antisymmetric. 
Explicitly, 

3 / 1 \ 3/4 . , 



$ N (ri,r 2 ,f 3 ;^) = ]J ( ^2 ) e " 2 * 2 ® [3]s=i/2,t=i/2® [l 3 ]c, (3.2) 

3 ( l \ 3/4 tf -m 2 

^A(ri,r 2 ,f 3 ;R)=Y[( — \ e <g> [3] s=3/2 , T=3/2 <g> [l 3 ] c , (3.3) 

n—1 ^ ' 




$N*(1440)(n, 7=2,7=3 = \ -01 - J -<p 2 ® [3] S= l/2,T=l/2 ® [1 ]c , (3.4) 



where [3]st stands for a completely symmetric ST/ (4) spin-isospin wave function and 
[l 3 ]c for a completely antisymmetric SU(3) colour wave function [63], and 

3 _ 2 



Mri,f 2 ,f 3 ;R) = x(^) E 



k=l 



3 (r fe - R) 2 



2 b 2 



2 f 1 \ ( — 
^(ri,f 2 ,f3 ;j R) = -- (^-tt) ^ (^-^.(rl-^ne-^. (3.6) 

3 Vtt4&2 y j<k=i , =1 



Once single baryon wave functions have been constructed, one can proceed to 
study the two-baryon wave functions. Again the symmetrization postulate forces the 
wave function of a system of n identical quarks to be totally antisymmetric under 
the exchange of any two of them. As single baryon wave functions are already 
antisymmetric, there is an important simplification to construct two-baryon wave 
functions: one needs the antisymmetrizer for a system of six identical particles but 
already clustered in two antisymmetric groups. The antisymmetrizer can then be 
written as [64] 

1 ( 3 6 \ 

■V), (3.7) 

+ This is not, for example, the case for the N*(1535) where the total spin of the particle is the result 
of coupling the intrinsic spin and the relative orbital angular momenta of the quarks. 
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where P^ is the operator that exchanges particles i and j, V exchanges the two 
clusters, and N is a normalization constant. P^ can be explicitly written as the 
product of permutation operators in colour (C), spin-isospin (ST) and spatial (O) 
spaces, 

Pij = I'! j P% T P?j • (3-8) 
Taking into account that any two-baryon state can be decomposed in a symmetric 
plus an antisymmetric part under the exchange of the baryon quantum numbers, one 
can write for a definite symmetry (specified by / even or odd) and projecting onto a 
partial wave to make clear the effect of the exchange operator [65,66]: 



LST 
B\B-2 



(R) 




LST 



(3.9) 



where S, T and L correspond to the total spin, isospin and orbital angular momentum 
of the two-baryon system. A is the six-quark antisymmetrizer described above. The 
action of V, appearing in the antisymmetrizer, on a state with definite quantum 
numbers LST is given by, 



V 



4>; 




LST 



J LST 



(3.10) 



LST 



Then, due to the (1 — V) factor in equation (|3.7|l . the wave function vanishes unless: 

L + S 1 + S 2 -S + T 1 +T 2 -T + f = odd. (3.11) 

For non-identical baryons this relation indicates the symmetry / associated to a given 
set of values LST . The non-possible symmetries correspond to forbidden states. For 
identical baryons, B\ = E>2, such as nucleons (note that / has to be even in order 
to have a non- vanishing wave function), one recovers the well-known selection rule 
L + S + T = odd. 

Certainly, the effect of quark substructure goes beyond the (1 — V) factor- 
appearing in the antisymmetrizer and it also appears through the quark permutation 
operator , whose effect can be analyzed in part in a simple way through the norm of 
the two-baryon system. This is a measure of the overlapping between the two baryons 
and it shows out the consequences of the Pauli principle. The norm of a two-baryon 
system B\B 2 is defined as, 

Making use of the wave function (|3.9|l one obtains: 

^b'bI(R) = ^l(R) - C(S, T, /; B Y B 2 ) N?{R) , (3.13) 



(3.12) 
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Table 2. Spin, isospin coefficients. 


"+" ("— ") refers to / even (odd). 


(S, T, f) C(S, T, /; NN) 


C(S, T, /; NA) 


C{S,T, /;AA) 


(0,0,+) 7/9 




1/3 


(0,1,+), (1,0,+) -1/27 




1/9 


(0,2,+), (2,0,+) 




-1/3 


(0,3,+), (3,0,+) 




-1 


(1,1,+) 31/81 


31/27 


1/27 


(1,1,-) 


1 




(1,2,+), (2,1,+) 


1/9 


-1/9 


(1,2,-), (2,1,-) 


-1/3 




(1,3,+), (3,1,+) 




-1/3 


(2,2,+) 


7/3 


1/3 


(2,2,-) 


1 




(2,3,+), (3, 2,+) 




1 


(3,3,+) 




3 



where J\f^(R) and A/"£ X (P) refer to the direct and exchange kernels, respectively. The 
direct kernel corresponds to the identity operator appearing in the antisymmetrizer, 
while the exchange kernel arises from the Py operator. C(S,T, /; B\B 2 ) is a spin- 
isospin coefficient defined as follows, 

C{S,TJ;B 1 B 2 ) = [(S 1 (123)S 2 (456)|P i f T |S 1 (123)S 2 (456))sT 

1 + Qb 1 b 2 

+ (-l)/( J B 1 (123)P 2 (456)|P l f T |P 2 (123)P 1 (456)) S T] . (3.14) 

These spin-isospin coefficients, summarized in table |21 determine the degree of Pauli 
attraction or repulsion as we will see later on. Let us note that for the NN*(1440) 
system, being the N* (1440) spin-isospin wave function completely symmetric, the spin- 
isospin coefficients are the same as for the NN case. The explicit expressions of the 
direct and exchange kernels depend on the baryons considered. For the NN, NA and 
AA cases (the NN*(1440) system is much more involved [67] and will be discussed 
below) the norm kernels are given by, 

3P^ 
"4 b 2 
3P^ 
'4 b 2 

where il are the modified spherical Bessel functions. 

To examine the physical content of Af, it is convenient to take the limit of the 
distance between the baryons approaching zero, R — > 0, 



A/HP) = 4tt exp - 



Afl*(R) = 4tt exp 



1 L 



IL 



3P^ 
4 b 2 

P2 
46 2 



(3.15) 



1 



3P? 
4 b 2 



1 



l-3---(2L + l) 



P? 

4b 2 



[Z L -C{S,TJ-B X B 2 ) 



C(S,T,f;B 1 B 2 )] 



(3.16) 



2(2L + 3) 

Of significant interest are those cases where 

3 L = C(S,T,f;B 1 B 2 ), 

because it implies that the overlapping of the two-baryon wave function behaves as 
P 2L+4 instead of the centrifugal barrier behaviour R 2L , indicating that quark Pauli 



(3.17) 
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Figure 1. NN*(1440) overlapping for L = partial waves. 



blocking occurs, the available spin-isospin-colour degrees of freedom saturate and then 
some quarks are Pauli expelled to higher orbits. This suppression in the overlapping 
of the two-clusters, that may be a source of short-range repulsion, is not present at 
baryonic level in those cases where B\ ^ B2, because of the distinguishability of 
baryons. Looking at table |21 we differentiate the following cases: 

• In the NN system there are not Pauli blocked channels. 

• In the NA system some partial waves present quark Pauli blocking, those 
corresponding to (S, T) — (1, 1) and (2, 2), with orbital angular momentum L = 
(/ = odd). Pauli blocking will translate into a strong short-range repulsion that 
can be checked experimentally looking at nd elastic scattering [68] . We will return 
to this point in section 14.2.21 

• In the AA system, the spin-isospin coefficients fulfill equation (|3.17(1 for the cases 
(S,T) = (2,3) and (3,2) both with orbital angular momentum L = (/ = 
even) . It is also important to mention the existence of quark Pauli blocking for a 
channel with L ^ which is a characteristic feature of the AA interaction, this 
corresponds to {S, T) = (3, 3) with orbital angular momentum L — 1 (/ — even). 
In a group theory language, 



showing the absence of the lower energy spatial states ([6]o for (S, T) — (2,3) 
and [51]o for {S,T) = (3,3)) forbidden by the symmetrization postulate. 

As previously said, for the NN*(1440) system the calculation is much more 
involved. For the most interesting cases, those where there is no centrifugal barrier, 
and in the limit R — > one obtains [67] : 



Quark Pauli blocked channels would correspond to /=odd and C(S, T, /; NN)=1, or 
/=even and C(S, T, /; NN)=3/7. From the values given in table [21 it is clear that 
there are no Pauli blocked channels. However, looking at figure ^ one can see how in 
those S*-wave channels forbidden in the NN case, (S, T) = (0, 0), (1, 1), the overlapping 



*aa[ 5 S 2 (T = 3)] = [2 3 ] c [42] [51]sr = *aa[ 7 S 3 (T = 2)] 
*aa[ 7 -P 2>3 ,4(T = 3)] = [2 3 ] a [33] o [6] st 



(3.18) 




(3.19) 
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gets suppressed as compared to the allowed channels (S, T) = (1,0), (0, 1). This is a 
remnant of the near to identity similarity of N and N*(1440), and it will have an 
important influence when deriving the NN*(1440) interaction. 

3.2. NN short-range repulsion 

As said above, there is no quark Pauli blocking in the NN system. Therefore the short- 
range repulsion experimentally observed requires a different microscopic explanation. 
More than twenty years ago a NN short-range repulsion was derived from an interplay 
between quark antisymmetry and OGE quark dynamics [26,27,30,32,33]. The origin 
of the NN short-range repulsion was understood in a simple and intuitive manner in 
terms of the energy degeneracy induced by the color magnetic hamiltonian between 
the different spatial symmetries in the two-baryon wave function [69]. Later on, the 
need to incorporate chiral potentials (OPE and OSE) in the description of the NN 
interaction forced a revision of the role played by the OGE dynamics and brought 
forth a new understanding of the short-range repulsion at the microscopic level. 

To settle this process let us first revise the OGE based explanation to establish 
the notation and physical arguments. In a group theory language, a completely 
antisymmetric six-quark state asymptotically describing two free nucleons in relative 
5-wave is given [for spin-isospin (S, T) = (0, 1) or (1,0)] by [31], 

* 6 ,(fl->oo) = * {6} + y^j v]/ {42} , (3.20) 

where 

* {6} = [2 3 ] c [6] [33] 5 t , 

* {42} = [2 3 ] c [42]o ([33]st ~ [51]st) , (3.21) 

the subindex C, O or ST indicating the colour, spatial or spin-isospin representation 
respectively. This six-quark wave function expressed in terms of the baryon-baryon 
basis would contain two-nucleon, two-delta and two-baryon coloured-octet states. 
Concerning the spatial part, L = may be obtained not only from six quarks in 
0s states ([6]o), but also from four quarks in 0s states plus two quarks in excited p 
waves ([42]o). When the distance goes to zero the {42} configuration represents an 
excited state. In fact, if we assume all quarks in a harmonic oscillator potential of 
angular frequency u, their wave functions will be given by 

/ 1 \ 3 / 4 

- T 3 t (57s) " 2 ik Y ^ ■ <3 - 22 » 

v TO= — 1 v v ' 

where we have assumed equal probability for all the third angular momentum 
components of the I = 1 quark orbital excited state rjo p (r). The excitation energy is 
then given by 

A ho ee [E {42} - E {6} ] ho = 2hu^ , (3.23) 
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the [6]o being the lowest in energy. This situation may be changed due to a particular 
dynamics. Let us for example analyze the case of the interaction via the OGE of 
reference [29] , 



1 



o> s Xi ■ A, 



2tt 
3 rriq 



aj £(r%) 



(3.24) 



where the values of the parameters are taken to reproduce the A— N mass difference: 
a s = 1.39, m q = 300 MeV and b — 0.6 fm. One can estimate the contribution of this 
interaction to the energy of the ^{42} or ^{6} configurations. First of all, taken into 
account that 



E x > 

i<j 



A, 



A; 



(3.25) 



where n is the number of particles, it is clear that quite approximately (1/Vy has been 
substituted by an average value) the Coulomb-like term of Voge does not contribute 

to the {42} — {6} splitting since ^JV \ij vanishes for a colour singlet. Hence, 
only the colour magnetic part (the term depending on the 5 function) has to be 
evaluated. Following reference [69], the matrix element for the {6} configuration [for 
(5,T) = (l,0)]is, 

8ci D 



< *{6} I - E Ci J & ■ 



{6} 



>= 



6 



= 48.83 MeV . 



(3.26) 



where the radial function aj is easily identified from l|3.24[) and c\ is calculated trough 



Cl = C®-' 08 =< T] Qs (f i )T] 0s (f j ) I Cj. 



nos(fi)r]os(f f j) >= 



12V2^6 3 m2 



D 
16' 



(3.27) 



D standing for the A— N mass difference (D = 293 MeV). Similarly, for the {42} 
configuration 

55c 2 



< * 



{42} 



E c y ( Xt ■ a j)^ 



{42} 



> 



-235.63 MeV, (3.28) 



where c 2 is the average two-quark interaction strength for the [42] q state, 



(•■2 



1 r 
15 



6 C, 



Os.Op 



379 
540 



Cl , 



(3.29) 



6, 8 and 1 denoting the number of the corresponding pairs. Finally one obtains 
Aoge = [^{42} — ^{6}] OGE = -284.46 MeV. The different character of the colour 
magnetic Hamiltonian for both configurations makes this difference to compensate the 
harmonic oscillator energy difference, Aho = 240.36 MeV, the [42] o spatial symmetry 
becoming the lowest in energy and almost degenerate with the [6]o- If the two spatial 
symmetry states are energy degenerate then for two free nucleons (R — > 00 in Eq. 
H3.20f> ) the [42] o dominates. To make the physics clear let us for a moment argue as 
if we had only the spatial [42] o symmetry (our results are not modified if the [6]o 
symmetry is included). At zero distance between the two nucleons the spatial [42] q 
symmetry implies a 2hu excitation in an oscillator basis. If the two nucleons are 
moved apart, in each nucleon the three quarks are in a spatially symmetric state and 
thus the configuration corresponds to a Ohio excitation. The 2Tluj energy has therefore 
to be in the relative motion. This means that for an S'-wave the relative motion must 
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have a node and thus cannot be described by a Os but instead by a Is relative wave 
function. In the asymptotic part of the wave function such a zero implies a phase shift 
as the one given by a hard-core potential [70] . 

Similar results are obtained with other OGE potentials. For example, using 
the parameters of reference [32] one obtains —291.08 MeV for the energy difference 
between the {42} and {6} configurations, while the harmonic oscillator energy gap 
is in this case 296.68 MeV. Once again the colour magnetic interaction produces a 
mixing of the symmetry states and a short-range repulsion. 

This explanation had to be revised in chiral constituent quark models where, in 
addition to the OGE, there are pseudoscalar and scalar Goldstone-boson exchanges 
between quarks. As a consequence the value of a s , which drives the OGE energy 
gap between the {42} and {6} configurations, is significantly reduced (due to the 
pseudoscalar contribution to the A— N mass difference) and correspondingly an energy 
degeneracy from the OGE is not attained. 

As a matter of fact if we recalculate the contribution of the OGE with the 
parameters of tabled we obtain 

< ^{6} I - c v & ' ^J'X^i ' S i) I *{6> >= 24 - 21 MeV > 

i<3 

< *{42} I - Ci J $ ' ^iX^< ' I ^i 42 } >= - 116 ' 78 MeV - ( 3 ' 3 °) 

Again, the colour magnetic interaction reduces the energy difference between the {42} 
and {6} configurations Aoge = —140.99 MeV, but the reduction of the energy gap 
is much smaller than the harmonic oscillator energy difference Aho = 309.10 MeV 
(note that the value of b in the CCQM is different than in pure OGE models). To go 
further let us analyze the OPE and OSE contributions. 

For the OPE potential l|2.5l) . the calculation of the corresponding spin-isospin 
matrix elements for a SU(4) irreducible representation is done in references [71,72]. 
The contribution of the OPE to the {6} configuration is, 

22 

< *{6} I X>i(*< • Zj-m ■ Tj) | * {6} >= lldi = — — ,—f. . , (3.31) 

where d\ — d^' * and the radial function dij can be identified from (|2.5p . For the 
{42} component we have, 

< [42] [33] S t | d » & ■ 5 iWi ■ f i) I t 42 lo >= 11 d, 

E 



< [42]o [51]st I 2^ d v ■ ■ f j) I I 42 ]o [51]st > = 35 d, (3.32) 

where 



15 



6d° s ' 0s + 8d° s ' 0p + d°f 0p 



(3.33) 



Then, from equation (|3.21() 

< *{42} I E d v & ■ ■ Tj) I *{42} > = 23 = -335.34 MeV. (3.34) 

i<j 

It is clear from the above expressions that differently than in the OGE case, the OPE 
contribution has the same sign for both configurations. Regarding the energy gap 
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between the spatial symmetries, one gets Aqpe = [-^{42} — -^{6}] OPE = —68.62 MeV. 
For the scalar potential l|2.6[) we get Aose = [^{42} — ^{6}] OSE = 133.11 MeV, its 
effect being just the opposite to the OGE one. 

Therefore, in chiral constituent quark models three effects conspire against the 
energy degeneracy of the spatial symmetry states. First, the small value of a s which 
lowers the contribution of the OGE. Second, the partial cancellation between the 
OPE contributions to the {42} and {6} configurations, and third the cancelling effect 
of the OPE+OSE potential with respect to the OGE. Putting all together one obtains 
Aqge + Aqpe + Aqse = —76.56 MeV which, differently than in the OGE models of 
references [29] and [32], is much smaller than the harmonic oscillator energy difference, 
A ho = 309.10 MeV. 

Although the perturbative separate one-channel calculation carried out should be 
only considered as a valuable hint (see for instance Ref. [33] for configuration mixing 
effects) the results obtained make plausible to conclude that in chiral constituent quark 
models there is not enough energy degeneracy to account for the NN short-range strong 
repulsion as a node produced by the [42] o spatial symmetry. 

Certainly this result depends on the value of the regularization cut-off mass A, 
that as has been explained controls the pion/gluon rate. The dependence on the 
cut-off mass A of the energy gap [-£{42} — -^{6}] generated by the OPE and OSE is 
presented in table One finds a small dependence of the results on small variations 
of A around reasonable values. It is worth to point out that the strong correlation 
among all the parameters does not allow for independent variations of them [73] . The 
strong coupling constant a s has to be recalculated for each value of A to reproduce 
the A— N mass difference. The new contribution of the OGE is also given in 

table El 

Since explicit NN calculations have shown that chiral constituent quark models 
have enough short-range repulsion to reproduce the experimental data [41,43,74] the 
question that immediately arises is where does the short-range repulsion come from?. 
To look for the origin of this repulsive character of the interaction one should go 
beyond the energy difference and calculate the specific contribution of the interaction 
for each symmetry. In order to see the repulsive or attractive character of each term 
of the potential in both spatial configurations one has to subtract twice (one for each 
nucleon) the corresponding nucleon self-energy, given by 

£oge = (N\V OG e\N) = -8 Cl = -72.63 MeV, 

£ope = (A|Vope|A) = 15di = -311.40 MeV, (3.35) 

£ose = (N\V OS e\N) = 3ei = -66.90 MeV, 

where e\ — e" 5 ' 05 , and is easily identified from Eq. (|2.(i|l . One then obtains for the 



Table 3. Variation of the energy difference in MeV between the {42} and {6} 
configurations. The harmonic oscillator energy gap is A ho = 309.10 MeV. 



A (fm 1 ) Aqse Aqpe Aqge Aqge + Aqpe + Aose 



3.7 


f23.72 


-64.21 


-152.40 


-92.89 


4.2 


133.11 


-68.62 


-140.99 


-76.50 


4.7 


140.98 


-71.98 


-131.13 


-62.13 


5.2 


147.49 


-74.77 


-122.39 


-49.67 


5.7 


152.98 


-77.07 


-115.98 


-40.07 


6.2 


157.65 


-79.07 


-110.15 


-31.57 
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Figure 2. Contribution of the different pieces of the interaction to the energies 
of the [42]o and [6]o spatial symmetries. Energies are given at R=0. 



{6} configuration, 

e oge+ope+ose = ( 169 47 + 486.95 - 200.70) MeV = 455.72 MeV, (3.36) 
and for the {42}, 

E {42f +OPE+OSE = ( 28 ' 48 + 418 - 33 ~ 67 - 65 ) MeV = 379 - 16 MeV - ( 3 - 37 ) 

The chiral potential produces strong repulsion in both cases due mainly to the OPE. 
In figure |21 it can be seen the effect of the different terms of the interaction. The 
OGE, repulsive in both configurations, reduces the energy difference whereas the 
OSE, attractive in both configurations, increases it. The OPE produces a strong 
repulsion in both symmetries. The net effect is an energy difference of about the same 
value obtained in the harmonic oscillator but with an additional repulsion in both 
symmetries originated mainly by quark antisymmetry on the OPE. Therefore, in this 
type of models the NN S*-wave hard-core like behaviour should be mainly attributed 
to the strong repulsion in the [6]o spatial configuration (more precisely there is a 
compromise between repulsion and kinetic energy difference in both configurations 
giving rise to some configuration mixing). 

These simple images of the origin of the NN short-range repulsion are confirmed by 
means of explicit RGM or Born-Oppcnheimer calculations of the NN relative motion 
based on OGE [26,27,30,32,33] or OGE plus Goldstone boson exchanges [42,43,47,74]. 
Such calculations have been done considering explicitly the NN and AA components 
(explicit hidden color-hidden color states were considered in Rcfs. [32,33]) and the 
mixing induced between them by the different interactions. 
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4. The baryon-baryon potential 

It has become clear in the last years the major role played by baryonic resonances, 
in particular the low-lying nuclconic resonances A(1232) and N*(1440), in many 
electromagnetic and strong reactions that take place in nuclcons and nuclei. This 
justifies the current experimental effort to study nucleon resonances in several facilities: 
TJNAF with a specific experimental program of clcctroexcitation of resonances, 
WASA in Uppsala to study NN — > NN7T7T reactions, etc. In this context the 
knowledge of the interaction involving resonances derived in a consistent way is of 
great relevance. The interaction between the nucleon and a resonance has been 
usually written as a straightforward extension of some pieces of the NN potential 
modifying the coupling constants extracted from their decay widths. Though this 
procedure can be appropriate for the very long-range part of the interaction, it is 
under suspicion at least for the short-range part for which the detailed structure of 
baryons may determine to some extent the form of the interaction. It seems therefore 
convenient to proceed to a derivation, besides the NN potential which will serve to 
fix the quark potential parameters, of the NN — > NR (R : resonance), NR — > NR, 
and RR — > RR interactions based on the more elementary quark-quark interaction. 
The main comparative advantage of the quark treatment comes out from the fact 
that as all the basic interactions are at the quark level, the parameters of each vertex 
(coupling constants, cut-off masses,...) are independent of the baryon to which the 
quarks belong, what makes its generalization to any other non-strange baryonic system 
straightforward. The other way around, the comparison of its predictions to the 
experimental data available serves as a stringent test of the quark potential model. 

The derivation of the dynamics of a two-baryon system from the dynamics of its 
constituents is a tough six-body problem whose solution cannot be exactly obtained 
even for the non-relativistic case. This forces the use of approximate calculation 
methods. In the literature two methods have been mainly used to get baryonic 
interactions from the dynamics of the constituents: the resonating group method 
(RGM) and the Born-Oppenheimcr (BO) approximation. We resume their most 
relevant aspects. 

4--1- Calculation methods 

4- 1.1. Resonating group method potential The RGM [75], widely used in nuclear 
physics to study the nucleus- nucleus interaction, can be straightforwardly applied to 
study the baryon-baryon interaction in the quark model. It allows, once the Hilbert 
space for the six-body problem has been fixed, to treat the inter-cluster dynamics in 
an exact way. 

The formulation of the RGM for a system of two baryons, B\ and i?2, starts from 
the wave function of the six-quark system expressed in terms of the Jacobi coordinates 
of the baryons. Then, the spatial part factorizes in a product of two three-quark cluster 
wave functions and the relative motion wave function of the two clusters so that, 

* BlB2 = A MP) 4>bM Bi ) <PbMb 2 ) xITb 2 £c[2 3 ]], (4.1) 

where A is the six-quark antisymmctrizer, ip(P) is the relative motion wave function 
of the two clusters, <pBi(p$, B . ) is the internal momentum wave function of baryon Bi, 
and are the Jacobi coordinates of the baryon Bi. xffs 2 denotes the spin-isospin 
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wave function of the two-baryon system coupled to spin S and isospin T, and, finally, 
£c[2 3 ] is the product of two colour singlets. 

The dynamics of the system is governed by the Schrodinger equation: 

{H - £t)|* >= => < 6V\(H -E T )\V >=0, (4.2) 

where 

6 ^2 6 
i=l q i<j = l 

Tqm being the center of mass kinetic energy, V qq the quark-quark interaction, pi 
the trimomentum of quark i, and m q the constituent quark mass. In equation 
(|4.2|) the variations are performed on the unknown relative motion wave function 
<p(P)- Assuming harmonic oscillator wave functions for ef>g i equation Q4.2[l . after the 
integration of the internal degrees of freedom of both clusters, can be written in the 
following way [76] , 




V D (P, Pi) + W EX (P, PM <p{P)dPi = , (4.4) 



where E = Et — Eg 1 — Eb 2 = Et — Ei n , Ei n being the internal energy of the two-body 
system and Vd(P, Pi) and W-ex(P,P) are the direct and exchange RGM kernels, 
respectively. Vb contains the effect of the interaction between baryonic clusters 
whereas Wex gives account of the quark exchanges between clusters coming from 
the identity of quarks. They are evaluated in detail in reference [77]. 

Note that if we do not mind how Vd and Wex were derived microscopically, 
equation 14. 4f) can be regarded as a general single channel equation of motion including 
an energy-dependent non-local potential given by the sum of Vd and Wex- This is 
the RGM baryon-baryon potential. 

4-1-2. Born-Oppenheimer potential The BO method, also known as adiabatic 
approximation, has been frequently employed for the study of the nuclear force from 
the microscopic degrees of freedom [25,78]. It is based on the assumption that quarks 
move inside the clusters much faster than the clusters themselves. Then one can 
integrate out the fast degrees of freedom assuming a fixed position for the center of 
each cluster obtaining in this way a local potential depending on the distance between 
the centers of mass of the clusters. The potential is defined in the following way [79], 

VB l B- 2 (LST)^BiBi{L> S'T){R) = Vf / T T (i?) - Vf / T T (oo) , (4.5) 

where 



with ^b S b T (R) given by equation (|3.9|l . The subtraction of V£ g T J '(oo) assures that 
no internal cluster energies enter in the baryon-baryon interacting potential. 
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Figure 3. Different diagrams contributing to the NN interaction. 



4-2. Results 

Both methods permit to evaluate the influence of the Pauli principle at the quark 
level on properties of the baryon-baryon interaction. The main conceptual difference 
between the resulting interactions is that the BO potential is local while the RGM 
one is non-local, 



t/RGM 
V B t B 2 - 

V B 1 B 2 - 



B3B4 



'B^Ba 



V B t B 2 - 

yBO 
V B,B 2 - 



B 3 B 4 C#) R') 

b s bA r )- 



(4.7) 



This means that the T matrix calculated solving a Lippmann-Schwinger equation 
has a different off-shell behaviour * and thus will give different results, the larger the 
difference the more the particles of the system under consideration explore the off-shell 
region. 

In both cases the calculation of a baryonic potential from the quark dynamics 
involves, due to the antisymmetry operator, the calculation of many different diagrams 
that are depicted in figure [3] for the case of the NN interaction. The direct terms are 
represented by diagrams (a) and (b), whereas quark-exchange terms correspond to 
diagrams (c) — (g). The number in square brackets corresponds to the number of 
equivalent diagrams that can be constructed (it counts the number of quark pairs 
which are equivalent to the pair ij singled out in the figure). Diagram (a) cancels 
almost exactly with the self-energy term Vf , s " s T T (oo), it gives a small contribution at 
short-range. Diagram (b) generates the asymptotic behaviour of the interaction. The 
relevance of the quark-exchange terms, diagrams (c)— (g), depends on the overlap of 
the baryon wave functions. They are responsible for the short-range structure of the 
quark-model-based potential and they vanish when the two baryons do not overlap. 
Next we present results for the NN, NA, AA, and NN*(1440) systems. 



* The on-shcll behaviour is very similar, in fact one can almost achieve on-shell equivalence by fine 
tuning the quark model parameters [80]. 
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Figure 4. NN OSE potential for two different spin-isospin channels. The dashed 
line represents the parametrization used in reference [41] . The solid line represents 
the chiral constituent quark model result [43]. 

4- 2.1. The NN interaction Given the huge amount of experimental data available 
on the NN system, a detailed calculation of the scattering phase shifts and bound 
state properties will be presented in section [SJ We use this section to discuss 
qualitative important aspects of the NN interaction that are naturally described in 
chiral constituent quark models. Let us first mention that the identity of quarks gives 
rise to the well-known selection rule L+S+T=odd (see section f5TTf) 

By construction the chiral quark pion potential reproduces the NN long-range 
interaction (section 12.2(1 . It has been already discussed in section 13.21 how the quark 
substructure of the nucleon allows to understand the short-range behaviour of the 

5- wave NN interaction. Concerning the medium-range attraction, figure 0] shows the 
scalar potential obtained in the chiral constituent quark model [43] as compared to 
a parametrization at baryonic level used in reference [41] to fit the NN experimental 
data. While at baryonic level different coupling constants are used: <? ct nn 

/4tt = 3.7 

for (S,T) = (0,1) and ^ nn /4tt = 1.9 for (S,T) = (1,1), the quark model result is 
obtained in both cases from the same chiral quark coupling constant. Moreover, the 
deuteron binding energy is also reproduced [62], while the scalar coupling constant 
used at baryonic level in this case is once more a different one, ff^NN/^ 71 " = 2-55 [81]. 

Another important feature of chiral constituent quark models is that although 
asymptotically the spin-isospin structure of the different terms of the baryon-baryon 
potential is the same as the corresponding terms of the quark-quark one, at short 
distances quark exchange generates a rather involved spin-isospin structure. This is 
due to the antisymmetrization operator, that can be factorized as A = (1 — 9P 3 6)(1 — 
V). When the quark spin-isospin operators are transformed algebraically to find the 
corresponding baryonic operators different structures are generated [64]. Taking into 
account that P 3 ^ T = j(l + c? 3 • <? 6 )(1 + f 3 • f 6 ), the effective NN interaction can be 
decomposed as follows, 

*nn T) = MR) + Vi{R) a Bl ■ a B2 + V 2 (R) f Bl ■ tb 2 + V 3 (R) a Bl ■ &b 2 f Bl • tb 2 , (4.8) 

where Vi(R) are functions of the interbaryon distance R, and a Bi {r Bi ) are the spin 
(isospin) baryonic operators. The Vi(R) functions can be obtained from the calculated 
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Figure 5. Spin-isospin independent part (Vb) and spin-isospin dependent parts 
(Vi, Vi, and V 3 ) of the OSE potential. 



potential for the different (S,T) channels, 
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(0,1 
NN 



(4.9) 



'NN 



The decomposition (14.81) comes out from any spin-isospin structure of the quark- 
quark potential. In figure [5] the spin-isospin independent, Vb, and the spin-isospin 
dependent, Vi, Vi and V3, terms generated by the OSE potential are shown. The 
main component of the interaction is, as expected, scalar and attractive, however a 
small spin-isospin dependence appears. This dependence is a completely new feature 
with respect to the usual scalar exchange at baryon level and may play a significant 
role in the understanding of different reactions such as the p(a, a') or p(d, d') [82]. 

The same decomposition can be applied to the OPE. In this case, the relative 
strength of the spin-isospin dependent and the tensor terms at the baryonic level is 
largely reduced with respect to the quark level case, equation (|2.5|) . This allows for a 
simultaneous explanation of the p(p, A ++ )n and p(n,p)n reactions [58] what is difficult 
to get when using baryonic meson-exchange potentials [83]. 



4-2.2. The NA interaction The inclusion of NA intermediate states has been 
considered as a possible improvement of NN interaction models at intermediate 
energies for a long time [84-86]. The NA interaction has been usually described by 
means of baryonic meson-exchange models [87] or parametrized by phenomenological 
potentials [68,88,89] with coupling constants and cut-off masses not well determined 
due to the lack of sufficient experimental information about the NA system. 

Since the A is not a stable particle but rather a 7rN resonance, one must establish 
what it is understood by the NA interaction. The A is considered as an elementary 
particle, being the coupling to the 7rN continuum the responsible for its width [90]. 
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Figure 6. 5 S 2 {T = 1) and 3 S 1 (T = 2) NA potentials. 

The vicinity of a nucleon modifies the properties of the A, because they can exchange 
particles between each other, eventually a virtual boson. This exchange, which looks 
like the interaction between two stable particles is what is considered as the NA 
interaction. The modification by the coupling to the continuum (A width) should be 
implemented when treating any particular problem. 

At the quark level, the NA interaction has been derived in the chiral constituent 
quark model [66] and used to study the NN system above the pion threshold [91]. 
Figure shows the NA potential calculated for two partial waves of different isospin, 
5 S , 2 (T = 1) and 3 Si (T = 2). In one case the contribution of the different terms of 
the potential has been separated. As in the other case this separation is qualitatively 
similar, the contribution of the different diagrams of figure is presented. The effect 
of quark antisymmetrization can be extracted by comparing the total potential with 
the term V36 [diagram (b) in figure y], the only significant one that does not include 
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Figure 7. 5 S 2 (T = 1) and 3 Si(T = 1) NA potentials. The solid line represents 
the quark-model based result, the dashed-dotted line the contribution of the OGE 
and the dotted line that of the OPE+OSE (indistinguishable of the total potential 
for the 3 5i(T = 1) partial wave). The dashed line corresponds to the baryonic 
meson-exchange potential of reference [87]. 
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Figure 8. a Si(T = 1) and 5 S2(T = 1) NA phase shifts and separable potentials 
of reference [68]. 

quark exchanges. All the exchange diagrams do not appreciably contribute beyond 1.5 
fm, where the overlap of the nucleon and A wave functions is negligible. Above this 
distance the interaction is driven by the V36 term and it equals the total interaction. In 
general, we see how the V12 -P36 term [diagram (c)] generates additional attraction and 
it is the V 36 P 36 term [diagram (d)] the main responsible for the short-range repulsion. 
The behaviour of the other diagrams depends on the partial wave considered. 

In section IXT1 the presence of NA quark Pauli blocked channels, specifically the 
3 Si(T = 1) and 5 S2{T — 2) was discussed. In figure we can see how quark Pauli 
blocking, existing in the S%(T = 1) partial wave, translates into a strong short-range 
repulsion, while this behaviour is not observed for non-Pauli blocked partial waves 
as it is the case for the 5 S*2(T = 1) channel. The strong repulsion is due to the 
combined effect of the fast decrease of the norm of the six-quark wave function when 
R — ► with the presence of non-vanishing direct terms in the potential. Actually, 
the repulsion manifests itself through the direct contribution of the OPE and OSE 
potentials. Let us note that if there were no direct contributions, i.e., only quark- 
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Figure 9. 3 P (T = 1) and 5 D (T = 1) NA phase shifts. The solid line represents 
the quark model results and the dashed line those of the baryonic meson-exchange 
model of reference [87]. Experimental data are from reference [92]. 



exchange diagrams contributed then the resulting behaviour would be quite different 
since the quark-exchange terms go to zero with the same power of R as the norm. 

For the sake of comparison we have also plotted in figure [7| the result for the 
baryonic meson-exchange model of reference [87]: the strong short-range repulsion 
does not appear. The only possibility to simulate it would be the use of large cut-off 
masses, but this would produce instabilities in the short-range part of the interaction. 
Moreover, the baryonic model gives the same structure in both partial waves, the only 
difference between them being a spin-isospin factor. 

In figure |SJ a separable NA potential reproducing the phase shifts (dots in figure 
IHJ obtained from the experimental ird elastic cross section is shown [68]. One can 
infer a hard-core in the 3 Si(T = 1) partial wave from the rapidly varying phase 
shifts changing sign for a pion energy of 219 MeV and a smooth behaviour for the 
5 S , 2(7 1 = 1) partial wave. This behaviour is quite similar to the one predicted by 
the chiral constituent quark model (see figure 0| . When not considering the quark 
substructure a baryon-baryon hard-core has to be introduced by hand to reproduce 
the experimental data [88,89]. 

Nevertheless, potentials are not observable quantities and one should use them 
to calculate phase shifts to be compared to data, whose extraction is still a matter of 
controversy. There are a few NA channels which are susceptible of being parametrized 
in terms of phase shifts. In figure [!|] we plot the quark- model result for two uncoupled 
isospin one NA partial waves, 3 Po and 5 Do- Although the error bars are still big, 
quark-model results agree reasonably well with the experimental data, whereas the 
baryonic meson-exchange predictions, due to its unsmooth character at short distances 
[91], are far from the data in the 3 Po case. 

4-2.3. The A A interaction The effect of A A components has been considered when 
studying NN phase shifts at intermediate energies and deuteron properties at the 
baryon level [93,94]. This treatment suffers from the same shortcomings mentioned 
in the NA case. 

At the quark level the two identical baryon selection rule, L-r-S+T=odd, comes 
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Figure 10. 3 Si(T = 2) and 3 Pi(T = 3) A A potential. 



out from quark antisymmetrization. From the chiral constituent quark model potential 
the AA interaction has been derived in reference [71]. The potential obtained is drawn 
in figure ^1] for two partial waves of different isospin, 3 Si(T — 2) and 3 Pi(T = 3). 
In one of them the contribution of the different exchanges has been separated and in 
the other the contribution of the dominant diagrams showed in figure EH are depicted. 
As can be seen, the dominant terms are diagrams (b) and (d), the others giving 
almost a negligible contribution. In both cases the direct contribution, diagram (b), 
(only generated by quark-meson exchanges) is attractive, while the effect of the quark- 
exchange diagram (d) (due to quark antisymmetry) is to generate repulsion. Beyond 
1.5 fm, the interaction is driven by the V36 term that equals the total interaction. 

According to the discussion in section l3~Tl the 5 S2{T = 3), 7 S^(T = 2), and 
7 P2,3,i{T = 3) channels correspond to quark Pauli blocked states. We compare in 
figure ITU the AA potential for two S waves, a Pauli blocked channel, 5 5 I 2(T = 3), and 
a non-Pauli blocked one, 5 S2(T = 1). While the last one presents a soft attractive 
behaviour at short distances, the 5 S , 2(T = 3) channel shows a strong repulsive core 
coming from the fast decreasing of the overlapping of the two-baryon wave function 
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Figure 11. 5 5 2 (T = 3) and 5 S 2 (T = 1) AA potential. 
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Table 4. Character of the AA interaction obtained in references [79,95], reference 
[98] , and reference [71] . The (*) denotes those channels where a bound state would 
be favoured. 



(S,T) 


Refs. [79,95] 


Ref. [98] 


Ref. [71] 


(0,1) 


Attractive 




Attractive (*) 


(0,3) 


Weakly repulsive 


Attractive 


Weakly repulsive 


(1,0) 


Attractive 


Repulsive 


Attractive (*) 


(2,3) 


Repulsive 




Repulsive 


(3,0) 


Attractive (*) 


Attractive (*) 


Weakly attractive 


(3,2) 


Repulsive 




Repulsive 



(denominator in l|4.6[l ) together with non- vanishing direct terms in the numerator. 

It is interesting to analyze the possible existence of bound states with the chiral 
constituent quark model since this has been the object of study of the AA system 
with quark models containing confinement plus OGE potentials [79,95-98]. In tabled 
results of reference [71], based on the chiral constituent quark model, are compared to 
those of references [79,95,98], based on a OGE model. With respect to the possibility of 
having a favoured bound state in the (5, T) — (3, 0) channel, as predicted in [79,95,98], 
reference [71] gets no OGE attraction but repulsion, although as it is shown in the table 
the total interaction is weakly attractive mainly due to the scalar exchange potential. 
Reference [71] predicts the (S, T) = (0, 1) and (1, 0) channels to be the most attractive 
ones. This attraction, combined with the fact that they are S waves and therefore 
the centrifugal barrier is not active, makes them the best candidates for possible AA 
bound states [99], as will be discussed in section 17. 2. til 

4.2.4. The NN*(1440) system The N*(1440) (Roper) couples strongly (60-70%) 
to the 7rN channel and significantly (5—10%) to the oN channel [100]. Its role in 
nuclear dynamics as an intermediate state has been analyzed at the baryon level. 
The presence of NN*(1440) configurations in the deuteron was suggested long ago 
[101-104]. Graphs involving the excitation of N* (1440) appear also in other systems, as 
for example the neutral pion production in proton-proton reactions [105] or the three- 
nucleon interaction mediated by 7r and tr exchange contributing to the triton binding 
energy [106]. The excitation of the Roper resonance has also been used to explain the 
missing energy spectra in the p(a, a') reaction [107] or the np — ► d(Tnr)° reaction [108]. 
Pion electro- and photoproduction from nucleons may take place through the N*(1440) 
excitation as well [109]. 

At the quark level the involved N*(1440) radial structure increases the number 
of diagrams contributing to the interaction. There appear diagrams generated by 
the two parts of the wave function, §\ and $2 in equation (|3.4() . Although involving 
interactions between excited and non-excited quarks, they can be classified as in figure 
121 The full set of diagrams for the NN*(1440) potential are explicitly given in [67] and 
those for the NN -> NN*(1440) potential in [110]. 

4-2.5. The NN* (1440) interaction In this case we shall classify the channels, for the 
reason that will be apparent in what follows, as forbidden channels, i.e., allowed in the 
NN*(1440) system but forbidden in the NN case, and allowed channels, i.e., allowed in 
both NN*(1440) and NN systems according to (|3.11|) . In figure IT21 the potential for a 
forbidden channel, 1 So (T=0), and for an allowed one, 1 Sq (T=l), are compared. As 
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Figure 12. Contribution of the different terms of the interaction to the 
1 S (T = 0) and 1 S (T = 1) NN*(1440) potentials. 



can be seen, forbidden channels in S waves are much more repulsive than allowed ones, 
this repulsion being mainly driven by the non-vanishing direct terms in the potential. 
Moreover, as detailed in reference [67], the potential for the forbidden 1 Sq(T = 0) 
channel is very much the same than for the allowed 1 P\(T = 0) and similarly for 
3 <Si(T = 1) and 3 Pj(T = 1) (in this last case with small dependence on J due to the 
tensor interaction) . This can be understood in terms of the Pauli and the centrifugal 
barrier repulsions. The Pauli correlations and the centrifugal barrier in the P waves 
prevent all the quarks to be in the same spatial state, much the same effect one has 
due to Pauli correlations in the S forbidden waves added to the presence of the radially 
excited quark in the N*(1440). 

The lack of available data and the absence of alternative quark model calculations 
to compare with makes convenient in this case to extract a pure baryonic potential 
whose difference with the total one emphasizes the effects of the quark substructure. 
The dynamical effect of quark antisymmetrization can be estimated by comparing 
the total potential with the one arising from diagram (b) in figure |21 V36 , which is 
the only significant one that does not include quark exchanges. The V36 potential 
turns out to be attractive everywhere. Let us note however that Pauli correlations are 
still present in the V36 potential, through the norm, in the denominator of equation 
(|4.6|l . To eliminate the whole effect of quark antisymmetrization one should eliminate 
quark-Pauli correlations from the norm as well. By proceeding in this way one gets 
a genuine baryonic potential, that will be called direct potential. The comparison 
of the total and direct potentials reflects the quark antisymmetrization effect beyond 
the one-baryon structure, see figure IT31 As V36, the direct potential is attractive 
everywhere. It becomes then clear that the repulsive character of the interaction for 
S and P waves at short distances is due to dynamical quark-exchange effects. For 
distances R > 2 fm the direct and total potentials are equal since then the overlap 
of the N and the N*(1440) wave functions is negligible and no exchange diagrams 
contribute appreciably. These results clearly illustrate that the use of a NN*(1440) 
potential as a generalization of the NN interaction should be taken with great care, 
specially for the forbidden channels. 

4-2.6. The NN — > NN* (1440) interaction The presence of two identical baryons 
in the initial state forces in this case the selection rule L+S+T=odd. For the 
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Figure 13. Comparison between the total and direct potentials for 1 5o(T = 0) 
and 1 S (T = 1) NN*(1440) potentials. 



NN — > NN*(1440) potential, most diagrams contributing to the interaction are due to 
the first term of the N*(1440) wave function J-JS, i.e., |[3](0s) 2 (ls)). In figure HI 
the potentials for So and 3 Si partial waves are shown (let us note that an arbitrary 
global phase between the N and N*(1440) wave functions has been chosen). Although 
the long-range part of the interaction (R > 4 fm ) comes dominated by the OPE, 
the asymptotic potential reverses sign with respect to both NN and NN*(1440) cases. 
Thus for S and D waves the NN — > NN*(1440) interaction is asymptotically repulsive. 
This sign reversal is a direct consequence of the presence of a node in the N*(1440) 
wave function what implies a change of sign with respect to the N wave function at 
long distances (for NN*(1440) there are two compensating changes of sign coming 
from the two Ropers). This is also corroborated by the study of the OSE interaction 
that is always asymptotically repulsive at difference to the NN and NN*(1440) cases. 
If the opposite sign for the N*(1440) wave function were chosen the long-range part of 
the interaction would be attractive but there would also be a change in the character 
of the short-range part. 

It is worth to remark that no quark antisymmetrization effects survive either in 
the numerator or in the denominator (norm) of equation l|4.6|) at these distances. In 
other words, the potential corresponds to a direct baryon-baryon interaction that can 
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Figure 14. 1 So and 3 Si NN -* NN*(1440) potentials. The long-dashed, dashed 
and dotted lines denote the OPE, OSE and OGE contributions, respectively. 
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be fitted as it is conventional in terms of a Yukawa function depending on the mass 
of the meson. 

The total potential turns out to be attractive from R = 1.5 fm down to a lower 
value of R different for each partial wave. This behaviour, related again to the node in 
the Roper wave function, contrasts with the elastic NN and NN*(1440) cases, where 
for instance for S and D waves the scalar part keeps always the same sign and gives the 
dominant contribution for R >0.8 fm. Below 0.6 fm the potential becomes repulsive 
in all partial waves. Nevertheless there are two distinctive features with respect to the 
elastic NN and NN*(1440) cases: in NN -> NN*(1440) the intensity of the repulsion at 
R =0 and the value of R at which the interaction becomes repulsive are significantly 
lower than in NN and NN*(1440) elastic potentials. This is a clear effect of the more 
similarity (higher overlap) in these cases between initial and final states what makes 
the Pauli principle more active. 

Let us also mention that at short distances, the interaction could be fitted in terms 
of two different Yukawa functions, one depending on the meson mass, m, the other 
with a shorter range depending on ^/(M N .( 1440 ) — Mn + m)m. These two Yukawa 
functions could be associated to the two diagrams with different intermediate states 
[toNN and mNN*(1440)] appearing in time ordered perturbation theory when an 
effective calculation at the baryonic level is carried out (let us realize that in a quark 
calculation the intermediate state is always mqq, the N*(1440)— N mass difference 
being taken into account through the N and N*(1440) wave functions). 



4.3. ttNN*(1440) and ctNN*(1440) coupling constants 

A main feature of the quark treatment is its universality in the sense that all baryon- 
baryon interactions are treated on an equal footing. This allows a microscopic 
understanding and connection of the different baryon-baryon interactions that is 
beyond the scope of any analysis based only on effective hadronic degrees of freedom. 
We will illustrate this discussion by means of the NN — > NN*(1440) transition 
potential, determining the 7rNN*(1440) and trNN*(1440) coupling constants. 

As has been discussed in section 14.2.11 asymptotically (R >4 fm) the OSE and 
OPE potentials have at the baryon level the same spin-isospin structure than at quark 
level. Hence one can try to parametrize the asymptotic central interactions as, 

t/OPE fm _l ffTrNN g7rNN*(1440) m v A 2 

'NN— >NN* (1440) v^V " 3 ^ ^ ^ 2(2Mr) A 2 _ ^ 

[{a N .a N )(f N .f N )} — — , (4.10) 
H 

and 

T/ OSE (m _ 9<?NN ff g NN*(1440) A 2 e~ m ' R 

^NN^NN* (1440) W - -J^ A 2 _ m 2 ~ R~ > < 4 1 ! 



where gi stands for the coupling constants at the baryon level and M r is the reduced 
mass of the NN*(1440) system. 

By comparing these baryonic potentials with the asymptotic behaviour of the 
OPE and OSE previously obtained one can extract the ttNN*(1440) and crNN*(1440) 
coupling constants in terms of the elementary irqq coupling constant and the one- 
baryon model dependent structure (see section I2.2JI . The sign obtained for the 
meson-NN*(1440) coupling constants and for their ratios to the meson-NN coupling 
constants is ambiguous since it comes determined by the arbitrarily chosen relative 
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sign between the N and N* (1440) wave functions. Only the ratios between 7rNN* (1440) 
and <rNN*(1440) would be free of this uncertainty. This is why only absolute values 
will be quoted, except for these cases where the sign comes as a prediction of the 
model. For this study the So will be used for simplicity. This is why only the central 
interaction has been written in equation (|4.10|) . 

The A 2 /(A 2 — m 2 ) vertex factor comes from the vertex form factor chosen at 
momentum space as a square root of monopole [A 2 / (A 2 +g 2 )] 1/ ' 2 , the same choice taken 
at the quark level, where chiral symmetry requires the same form for pion and sigma. 
A different choice for the form factor at the baryon level, regarding its functional form 
as well as the value of A, would give rise to a different vertex factor and eventually to 
a different functional form for the asymptotic behaviour. For instance, for a modified 
monopole form, [(A 2 — m 2 )/(A 2 — q 2 )] 1 ^ 2 , where m is the meson mass (m w or m CT ), the 
vertex factor would be 1, i.e. (A 2 — m 2 )/(A 2 — m 2 ), keeping the potential the same 
exponentially decreasing asymptotic form. Then it is clear that the extraction from 
any model of the meson-baryon-baryon coupling constants depends on this choice. We 
shall say they depend on the coupling scheme. 

For the OPE with A = 4.2 fin -1 , A 2 /(A 2 - m 2 ) = 1.03, pretty close to 1. As 
a consequence, in this case the use of this form factor or the modified monopole 
form at baryonic level makes little difference in the determination of the coupling 
constant. This fact is used when fixing g 2 qq /4ir from the experimental value of 
g 2 NN /47r extracted from NN data. To get S7rNN*(i440)/v47r we turn to the numerical 
result for the 1 S'o OPE potential and fit its asymptotic behaviour (in the range 
R : 5 — >9 fm, see figure IT5|l to equation H4.10JI . obtaining 



coupling constant is well defined. In reference [111] a different sign is obtained what is 
a direct consequence of the different global sign chosen for the N*(1440) wave function. 
The coupling scheme dependence can be explicitly eliminated comparing (77rNN*(i44o) 



i.e. 3 W NN*(i440)/v / 47r 
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Figure 15. Asymptotic behaviour of the OPE and OSE 1 So NN -> NN*(1440) 
potential (solid line), the dashed line represents the fitted curves according to 
14. 1 Ut and 14.111 . For the OSE potential both lines are indistinguishable. 
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with (/ttnn extracted from the NN potential within the same quark model, 

<?ttNN*(1440) 



0.25. (4.13) 

SttNN 

By proceeding in the same way for the OSE potential, i.e. by fitting the potential to 
equation (|4.11l) . and following an analogous procedure for the NN case one can write 

<7<tNN*(1440) 



= 0.47. (4.14) 

ffcrNN 

The relative phase chosen for the N*(1440) wave function with respect to the 
N wave function is not experimentally relevant in any two step process comprising 
N*(1440) production and its subsequent decay. However it will play a relevant role in 
those reactions where the same field (tt or a) couples simultaneously to both systems, 
NN and NN*(1440). In these cases the interference term between both diagrams would 
determine the magnitude of the cross section [107]. 

The ratio given in l|4.13|l is similar to that obtained in reference [111] and a factor 
1.5 smaller than the one obtained from the analysis of the partial decay width. Indeed 
one can find in the literature values for / w nn*(i44o) ranging between 0.27—0.47 coming 
from different experimental analyses with uncertainties associated to the fitting of 
parameters [108,109,112]. Regarding the ratio obtained in (|4.14(l . it agrees quite well 
with the only experimental available result, obtained in reference [107] from the fit of 
the cross section of the isoscalar Roper excitation in p(a, a') in the 10—15 GeV region, 
where a value of 0.48 is given. 

Furthermore, a very definitive prediction of the magnitude and sign of the ratio 
of the two ratios is obtained, 

g7rNN*(1440) _ q 53 gg-NN*(1440) ,^ ^ 

SttNN SctNN 

which is an exportable prediction of the chiral constituent quark model. 

For the sake of completeness we give the values of 5<jnn*(1440) an( i ffo-NN, though 
one should realize that the corresponding form factor A 2 /(A 2 — m 2 ) = 2.97 differs 
quite much from 1. Extracting the quark model factor dependence from the coupling 
constant (e m ^ 2 / 2 ) [67], that one may consider included in the baryon form factor, the 
results obtained 

5^nn* (1440) I — 1.14, that compares quite well with the value 
given in reference [107], 5^NN*(i440)/^ 7r = 1-33, and <7 2 NN /47r = 5.06. These coupling 
constants have been also determined in reference [113]. The results reported there are 
sensitive to the decay width of the sigma into two pions and the mass of the sigma 
as reflected in the large error bars given. Both quantities are highly undetermined 
in the Particle Data Book [100], the mass of the sigma being constrained between 
400-1200 MeV and its width between 600-1000 MeV. These values have been fixed 
in reference [113] to m a =500 MeV and r CT =250 MeV. Varying the mass of the 
sigma between 400 and 700 MeV for a fixed width of 250 MeV, the coupling constant 
according to equation (9) of reference [113] varies between 0.18—2.54. Taking a width 
of 450 MeV the resulting coupling is 0.27—1.64. In both cases, our values may be 
compatible with the N*(1440) decay and production phenomenology. 
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5. The NN system 

The most accurate description of NN scattering data and deuteron properties has 
been done in terms of boson exchanges at the baryon level in the framework of the 
old-fashioned time ordered perturbation theory, viz., the Bonn potential [114]. On a 
microscopic base different approaches have attempted to describe the NN interaction. 
Effective theories are very efficient to reproduce the high orbital angular momentum 
(L) partial waves but for low L they need to introduce a large number of parameters 
(more than twenty in reference [115]). Bag models found severe problems in describing 
the two-baryon dynamics due, on the one hand, to the critical difficulty of separating 
the center of mass motion and, on the other hand, to the unphysical sharp boundary 
of the bag and the problem of how to connect six-quark dynamics with external 
NN dynamics. In the so-called hybrid model approaches short-range dynamics based 
on the quark substructure combines with a long-range part given by either meson 
exchanges at the baryonic level (for example references [116] and [117] used the 
long-range interaction from the one- and two-pion exchange Paris potential [118]) 
or phenomenological potentials fitted to the experimental data [119]. Although they 
had a relative success in reproducing the phase shifts these hybrid models are not fully 
consistent. 

Chiral constituent quark models allow to incorporate the medium and long-range 
dynamics from the quark substructure in a natural way. They were first applied to 
the non-strange baryon-baryon interaction [42,43] and later on to derive the potential 
between all members of the baryon octet [47] . In the last case the model parameters for 
the exchange of the two pseudoscalars (pion and kaon) and the full scalar and vector 
octet mesons were independently fitted to the experimental data. This procedure, 
although very effective to fit the data, does explicitly break the chiral construction 
of the potential and masks on the model parameters some important physical effects, 
like the coupling to NA channels [120]. In the following we will refer to the model of 
reference [43], that has been detailed in sectional although we will use results of the 
quark model approach of Ref. [47] for comparison. 

5.1. Low-energy scattering parameters and deuteron properties 

The scattering length and the effective range provide a useful way to parametrize 
information on low-energy NN scattering because they are extremely sensitive to 
small variations on the strength of the NN force. Furthermore, these parameters 
can be related to other observables such as the deuteron binding energy. Therefore, 
the correct prediction of the scattering length, the effective range and the deuteron 
properties should be the first test for any model of the NN system to satisfy. Table 
resumes the results obtained by means of the chiral constituent quark model [80] 
compared to the quark model approach of Ref. [121] and some standard NN potentials 
[114, 122] and experimental data. For a correct description of the 1 Sq observables it 
is necessary to take into account the coupling to the 5 Do NA channel, which provides 
an isospin-dependent mechanism generating additional attraction. Hence a similitude 
between the 1 So and the 3 Si NN partial waves, both getting the necessary attraction 
by the tensor coupling to a I? wave ( 3 Di for 3 Si) [120] comes out. Besides it is worth 
to mention that the description of the low-energy scattering parameters requires a 
charge symmetry breaking term [114, 123] that has been taken into account by a slight 
modification of the chiral coupling constant. In the same table the results for the 
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Table 5. NN properties. a* a (t) is the singlet (triplet) scattering length and r s t t \ 
is the singlet (triplet) effective range. is the deuteron binding energy, Po 
the D— state probability, Qd the quadrupole moment, Ag the asymptotic 5-statc 
normalization, and n the D/S-state ratio. 

CCQM [80] Rcf. [121] Nijm II [122] Bonn B [114] Exp. 



Low-energy scattering parameters 





a s (fm) 


-23.759 


-23.76 


-23.739 


-23.750 


-23.74±0.02 




r s (fm) 


2.68 


2.58 


2.67 


2.71 


2.77±0.05 


3 s 1 


a t (fm) 


5.461 


5.399 


5.418 


5.424 


5.419±0.007 




n (fm) 


1.820 


1.730 


1.753 


1.761 


1.753±0.008 


Deuteron properties 


E d 


(MeV) 


-2.2242 


-2.225 


-2.2246 


-2.2246 


-2.224575 


Pd 


(%) 


4.85 


5.49 


5.64 


4.99 




Q d 


(fm 2 ) 


0.276 


0.270 


0.271 


0.278 


0.2859±0.0003 


As 


(fm- 1 / 2 ) 


0.891 




0.8845 


0.8860 


0.8846±0.0009 


V 




0.0257 


0.0253 


0.0252 


0.0264 


0.0256±0.0004 



deuteron properties are shown, presenting a good agreement with the experimental 
data. 

5.2. Deuteron configurations 

Although the deuteron has been usually described as a NN isospin singlet in even 
partial waves, i.e., 3 Si and 3 Di, it could also be understood as a linear combination of 
pairs of baryonic resonances provided they have the adequate quantum numbers. Such 
a description finds its natural framework in quark models for two reasons. First, from 
a quark model point of view baryon resonances are no more than internally excited 
nucleons. These excitations may occur even at low energies due to the virtual process 
NN — > NN*(N*N*) — > NN involving intermediate N*'s. Second, all baryon-baryon 
interactions can be generated from the same underlying quark-quark potential and, 
therefore, described in a consistent way. Furthermore, the presence of colour for the 
quarks adds the possibility of having new exotic components such as hidden colour 
state configurations: two colour octets coupled to a singlet. Although in some works 
their contribution to the deuteron wave function has been found to be as large as 
5% [124], it is possible to demonstrate that any hidden colour state can be expressed 
as a linear combination of physical states, i.e., baryon-baryon states [125]. 

The most prominent low-lying even-parity nucleon resonances are the A and 
the N*(1440). Being the deuteron isoscalar the NA configuration is forbidden 
and therefore AA and NN*(1440) components would be the relevant non-nuclconic 
configurations. For the AA configuration, reference [93] established a significant 
upper limit of about 0.4%. Moreover, although there are no experimental data on 
the NN*(1440) configuration, this has been advocated long ago to understand elastic 
proton-deuteron backward scattering at energies above pion threshold [126] or the 
angular distribution of deuteron photodisintegration at energies above E 1 = 100 
MeV [127]. Recent calculations have renewed the interest on these non-nucleonic 
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Table 6. Probability of the different deuteron components (%) 





NN 


NN*(1440) 




AA 






3 D 1 




:i D 1 


3 Sl 






95.3780 


4.6220 












95.1989 


4.5606 






0.1064 


0.0035 0.1243 


0.0063 


95.1885 


4.5377 


0.0022 


0.0148 


0.1224 


0.0036 0.1245 


0.0063 



components as they could be indirectly observed in several reactions as for example 
antiproton-deuteron annihilation [128], subthreshold antiproton production [129] or 
pd — > dp processes [130]. 

There are several deuteron multichannel calculations using different approxima- 
tions. Reference [131] proposed a formulation in terms of quark-shell configurations, 
later on projected onto physical channels. Reference [132] studied the effective num- 
bers for different resonance configurations making use of baryon wave functions ob- 
tained from the diagonalization of a quark-quark interaction containing gluon and 
pion exchanges in a harmonic oscillator basis including up to 2Hu> excitations, and 
deuteron wave functions obtained from the Paris potential. An upper limit of 1% for 
AA components and 0.1% for NN*(1440) was obtained. Possible effects of six-quark 
bags at short-distances were considered in [133]. In reference [134] the influence of 
N and A resonances on the NN interaction was studied. The most significant con- 
tribution was inferred from channels involving N and A ground states, although a 
quantitative calculation of the non-nucleonic configurations was not performed. 

The most complete multichannel calculation was done in [135] and [136] using 
a chiral constituent quark model. The different configurations and partial waves 
included in the calculation are shown in table El together with their probabilities 
(for instance the first row contains only the NN configurations 3 Si and 3 Di). In 
all cases the deuteron binding energy is correctly reproduced, being E4 = —2.2246 
MeV. The first remarkable result (the last row of table EJ) is that the probabilities of 
NN*(1440) channels are smaller than the AA ones being these later ones compatible 
with the experimental limits given in [93]. As compared to other values in the 
literature, the total probability of AA components is three times smaller than the 
one reported in [137], although presenting the same distribution with respect to the 
different channels. Regarding the calculations of reference [138], a similar distribution 
is obtained but the total probability being twice as much. Nonetheless, the comparison 
with these results should be done with care. For example, reference [137] carries 
out a coupled-channel calculation using for the NN — > AA transition potential a 
combination of n and p exchanges, and the Reid soft core potential for the rest. 
Such calculations are not fully consistent since they use NN potentials (Reid soft core 
or Paris) which were designed, in principle, without explicit A degrees of freedom. 
Therefore, parts of the contribution from the A degrees of freedom are implicitly 
included in the parameterization of the potentials. 

Regarding the NN*(1440) probabilities, they show a larger component of the 3 £>i 
partial wave against the 3 Si partial wave, in agreement with the ordering obtained 
by other calculations [104,134]. This can be understood taking into account that 
the tensor coupling, which is the main responsible for the presence of non-nucleonic 
components on the deuteron, is much stronger for the 3 5i(NN) — > 3 Di[NN*(1440)] 
transition than for the 3 L>i(NN) — > 3 5i[NN*(1440)] one, enhancing in this way the D- 
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wave influence with respect to the S'-wave component. Concerning the absolute value 
of the probabilities, they are a factor ten smaller than those reported in reference [104], 
where an estimation of 0.17% for the NN*(1440) configuration was obtained (0.06% 
for the 3 Si and 0.11 % for the 3 D\ partial wave). The difference with the result of 
the chiral constituent quark model may be understood in the following way. The 
deuteron is calculated using the NN Reid hard-core potential. When including 
NN*(1440) components, the channel coupling induces an attractive interaction on 
the NN system that needs to be subtracted out. Such a subtraction is done by 
reducing the intermediate range attraction of the central part of the Reid hard-core 
potential without modifying the tensor part (as done in reference [139] to calculate 
the probability of A A components). So, the strength of the tensor coupling to the 
NN*(1440) state is not strongly constrained. In fact, a variation of this strength may 
be compensated by a corresponding modification of the intermediate range attraction. 
Thus, the balance of these two sources of attraction cannot be disentangled in a 
clearcut way. This is a similar problem to the one arising in the NN partial wave 
when the coupling to the NA system is included. The same attractive effect may be 
obtained by a central potential or by a tensor coupling to a state with larger mass, 
being necessary other observables to discriminate between the two processes [140]. 

Although the probability of AA components is small, its influence for some 
specific observables can be important. Figure shows the momentum distribution 
k 2 r](k) (rj(k) = ^(fc)] 2 is the probability density in momentum space) for the 5*-wave 
deuteron components. For the 3 Si(NN) case, the solid line represents the results 
including all NN and AA partial waves while the dashed line refers to the calculation 
only with NN channels. For the 3 Si(AA) case the dotted line represents the result 
obtained with all NN channels and only the 3 <Si(AA) partial wave, while the dashed 
line includes all NN channels and the most important A A partial waves, 3 Si and 7 -Di, 
being almost identical to the full calculation. As can be observed, the AA components 
extend to higher momenta regions. This may influence the structure function B(q) 
which presents a zero for momentum around 7.1 fm _1 . Including only NN components, 
this result always comes lower, being the AA component a possible candidate to solve 
the problem [138,141]. 
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5.3. NN phase shifts 

A calculation of the NN phase shifts below the pion threshold is presented in this 
section. It includes coupling to AA channels for iso-singlet (T = 0) partial waves, 
and to NA and AA channels for iso-triplet (T = 1) partial waves. Experimental data 
have been obtained through the interactive program SAID [142] corresponding to the 
neutron-proton solution SP98. 

For the discussion of the results the phase shifts are divided into three groups. 
First, the L = partial waves. They are the most sensitive to the short-range part of 
the interaction and therefore one would expect that quark-exchange effects played for 
them an important role. A second group will be P, D and F waves. They are still 
sensitive to the short-range part and therefore to quark exchanges, but the middle- 
range and spin-orbit terms will also play a relevant role. Finally, those partial waves 
with orbital angular momentum L > 3 depend basically on the long-range part of the 
interaction. Although quark exchanges are not expected to be relevant for them, the 
correct description of these phase shifts supposes a crucial test to the chiral symmetry 
hypothesis used to derive the basic quark-quark interaction. 

5.3.1. S partial waves Figure IT71 shows the results for the and 3 «Si partial waves. 
A perfect fit is obtained for 1 Sq up to 250 MeV laboratory kinetic energy and up to 
300 MeV for 3 Si. As can be checked, the chiral components of the NN interaction 
(dashed line) do not provide enough attraction to reproduce the experimental data. 
As previously mentioned, in the channel the required attraction is supplied by 
addition of a coupling to the 5 D NA channel (dotted line). A complete agreement 
with the experimental data is obtained when the coupling to AA channels is included 
(solid line) . This result shows that the requirements of chiral symmetry (the presence 
and parameters of the scalar exchange between quarks) are fulfilled by the data 
provided that one includes the necessary additional physics into the problem. 

For the iso-triplet partial wave 3 Si, the coupling to A A channels has a smaller 
influence on the phase shift as can be seen on the figure, and can be simulated by a 
fine tuning of the mass of the scalar boson [135]. 

Let us emphasize that it is the strong spin-isospin independent repulsion generated 
by quark antisymmetry effects on the OPE what allows to reproduce the behaviour 
of the S'-wave phase shifts. 
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5.3.2. P, D and F partial waves Among the four L = 1 partial waves, only the Pi 
is not affected by the spin-orbit interaction. As seen in figure^|the chiral constituent 
quark model result (solid line) is in perfect agreement with the experimental data. 
It is interesting to notice that for this partial wave quark-exchange effects still play 
an important role. Indeed, in the same figure, the calculation done by removing the 
terms coming from quark antisymmetry is shown (dashed-dotted line), and the phase 
shift becomes attractive. This result is easily understood in terms of the direct OPE 
interaction (the dominant OGE interaction coming from the spin-spin term does not 
contribute to P- waves because of its S- like behaviour). At short-range the direct OPE 
potential is repulsive for S waves, but it is attractive for the 1 P\ wave due to the 
different sign of the spin-isospin matrix element. Only when quark antisymmetry 
terms are considered the OPE produces the correct experimental behaviour. 

In the same figure the 3 Pj triplet phase shifts are shown (notation as in figure lT7|) . 
One observes, as has been explained for S'-waves, how the coupling to AA channels in 
higher angular momentum iso-triplet partial waves is very small. The coupling to NA 
channels does not improve the description. The failure to reproduce the data may be 
associated to the lack of spin-orbit interaction. 

Some selected D waves are shown in figure IT§1 (notation as in figure IT7|l. As a 
general trend, one can say that the model gives too much attraction, except for the 3 Di 
partial wave where the experimental phase shifts are perfectly reproduced. While for 
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iso-triplet partial waves the coupling to NA channels is still important (worsening the 
quality of the results), the coupling to A A channels does not produce any considerable 
effect. For the 1 Z?2 partial wave the result obtained switching off the antisymmetry 
effects is again displayed by a dashed-dotted line, showing clearly how the influence 
of antisymmetry diminishes as compared to S and P waves. It is important to note 
that the D- waves are correctly reproduced up to about 80-100 MeV, which may be 
considered as the characteristic window where the NN interaction is basically governed 
by chiral symmetry [143]. Higher energies require an improvement in the description 
of short-range effects, as could be the spin-orbit force which still plays an important 
role for D waves. 

The same figure EH includes some selected F-w&ve phase shifts. In general there 
is a better agreement than in D-waves up to higher energies of 150-200 MeV. For 
iso-triplet partial waves the coupling to NA channels produces too much attraction, 
as it was the case for P and D waves, but the description is of the same quality as 
that of effective theories [143,144]. A comment about the 3 i<2 phase shift is in order. 
This wave is coupled to the 3 P2 wave, which has a strong influence of the spin-orbit 
interaction. The spin-orbit interaction is attractive for the 3 P2 wave and repulsive for 
the 3 i 7 2- Therefore, one expects that when the 3 P2 wave approaches the experimental 
data, the 3 P2 will also do the same. 
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Tkb < MeV > T |ab (MeV) 

Figure 20. J = 5 NN phase shifts. 



5.3.3. G, H and I partial waves In figure |2*U1 J = 5 G, H and / phase shifts are 
shown (notation as in figure ITTjl . The results obtained are in good agreement with 
the experimental data, being fundamentally dominated by the interaction of chiral 
origin. Iso-triplet partial waves show that the coupling to NA channels, so important 
for low orbital angular momenta, becomes smaller but contributes to improve the 
results. Once again they are similar to those obtained by effective theories [143, 144], 
indicating that these waves are governed exclusively by chiral symmetry. 

5.3.4- Mixing parameters Figure |^ shows the mixing parameters for the tensor 
coupling between different NN partial waves: 3 Si - 3 Di, 3 P 2 - 3 F 2 , 3 D 3 - 3 G 3 and 
3 G5 — 3 /5. The dashed line shows the results when only NN components are included. 
As can be seen, ei (which represents the tensor coupling 3 Si — 3 D\) and e 2 (which 
represents the tensor coupling 3 P 2 — 3 F 2 ) are not correctly reproduced. For e 2 one 
cannot draw any definite conclusion, because the model does not accurately describe 
the 3 P 2 and 3 F 2 phase shifts. However, despite the very good description of the 3 Si and 
3 D\ phase shifts, ei has the correct behaviour only at very low energies. This seems 
to indicate the need for a stronger short-range tensor force. Such interaction might 
be obtained from a simultaneous exchange of pions and gluons between quarks [145]. 
Higher angular momentum mixing parameters are governed by the OPE tensor term 
and present a good agreement with the experimental data. 



Quark-model study of few-baryon systems 



43 




100 200 300 100 200 

T lab (MeV) T |ab (MeV) 
Figure 21. Mixing parameters. 



5.3.5. Alternative quark model approaches Other constituent quark model 
approaches have been used to describe the NN dynamics. Ref. [74] used a generalized 
chiral symmetry restoration scheme giving rise to four chiral fields instead of the 
simple image of the 7r and a mesons. Similar results to the ones presented here are 
obtained. In Refs [47,48,121] Fujiwara and collaborators have performed a complete 
analysis of the NN scattering phase shifts and mixing coefficients. Their quark model 
approach incorporates a non-relativistic kinetic energy term, a quadratic confinement 
potential, the full Fermi-Breit interaction with explicit quark-mass dependence and 
the full octet of scalar and pseudoscalar meson-exchange potentials between quarks. A 
refined version of the model considered also vector meson-exchange potentials and the 
momentum-dependent Bryan-Scott term in the scalar and vector meson potentials. An 
impressive good agreement with the experimental data has been obtained for almost 
all partial waves and also for the mixing coefficients. However consistency with the 
one-body problem is lost (see also section 1775)) . The main difference with respect to 
the CCQM used here comes from the presence in the CCQM of a form factor as a 
consequence of the momentum dependent constituent quark mass generated by chiral 
symmetry breaking, Eq. (|2.3I) , that allows to regularize the non-central contributions 
of the OPE without resorting to multiple meson exchanges as in Refs. [47,48,121]. 
In fact, in these works an important contribution from the non-central terms of the 
octet meson-exchange potential is obtained. In particular, the quadratic spin-orbit 
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component cancels the strong one-pion tensor force. The interest reader could find a 
detailed description in the quoted works. 

5.4- The NN spin- orbit force 

The study of P and D phase shifts makes evident a lack of spin-orbit interaction in 
chiral constituent quark models. The origin of the spin-orbit force is a long-standing 
problem of the quark model in its description of the NN interaction and the baryon 
spectrum [19,41,52,64, 146-149]. On the one hand the spin-orbit force is known to 
be short ranged and thus a priori it may seem easy to understand on a quark model 
scheme, but on the other hand, it is a relativistic correction and therefore difficult 
a priori to accommodate in a non-relativistic approach. In baryonic meson-exchange 
models [114, 150] the spin-orbit force has its origin on the ui- meson exchange, also 
responsible for the short-range repulsion. In quark models the inclusion of vector 
meson exchanges between quarks belonging to different clusters provide contributions 
similar to the quark-exchange mechanism and would lead to double counting [57], as 
explained in section 12.11 Therefore, vector mesons are not included and one has to 
look for a different source of the spin-orbit force. 

The non-relativistic reduction of the one-gluon exchange diagram provides at 
(Pq/m q ) 2 order two different types of spin-orbit terms, a Galilei invariant (symmetric) 
term and a Galilei non-invariant (antisymmetric) term. The Galilei-invariant spin- 
orbit term has received considerable attention in the literature [41,147-149,151], 
because it presents the same structure as the standard spin-orbit interaction between 
nucleons. Using this source of spin-orbit force it was demonstrated in reference [151] 
that it is possible to reproduce the observed spin-orbit splitting in the NN interaction 
using a large quark-gluon coupling constant a s ~ 1.6. However, an effective meson- 
exchange potential between nucleons and not between quarks was used. When quark- 
meson exchange is incorporated instead, a s is considerably reduced and the resulting 
spin-orbit force is not large enough to reproduce the experimental P-wave splittings. 
In fact, when gluon and pion exchange are consistently treated at the quark level, the 
OGE spin-orbit force has to be multiplied by some numerical factor (usually between 
5—12 depending on the set of parameters [41,43,152]) in order to reproduce the 
experimental splittings in the NN phase shifts, with the counterpart that it severely 
disturbs the description of the excited negative-parity baryon spectrum [19, 151]. 
Regarding the antisymmetric spin-orbit OGE term, the strength of the full triplet- 
odd NN spin-orbit potential is greatly reduced and its sign may be even reversed 
when considered. This is the reason why it has not been usually considered in the 
literature, although there is no fundamental reason to justify this. 

Another source of spin-orbit force in chiral constituent quark models is the 
scalar exchange. In reference [73] it was shown that it is not enough to explain the 
experimentally observed effects. The strength of the NN one-sigma exchange spin-orbit 
force (OSE-LS) can be compared to the one-gluon exchange symmetric spin-orbit force 
(OGE-LS) in the adiabatic approximation. Figure 1221 shows the radial part of the NN 
adiabatic spin-orbit potential for 3 Pj and 3 Dj waves. The OGE-LS is generated only 
by quark-exchange diagrams what makes it highly non-local, while the OSE-LS has 
a non-vanishing direct term making it longer-ranged. The sign of the OSE-LS is the 
same for P and D waves, because the direct term dominates. However, the OGE-LS 
force has a different sign for P and D partial waves, since only the exchange terms 
contribute. 
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Figure 22. The adiabatic OGE and OSE spin-orbit potential of reference [73] 
for 3 Pj and 3 Dj partial waves. 



Finally, the spin-orbit potential generated by confinement is strongly dependent 
on the particular confinement model. For example, while the symmetric spin-orbit 
term associated with a scalar two-body confinement potential cancels the OGE 
symmetric spin-orbit force in the NN interaction [52], the one coming from a flip- 
flop confinement interaction interferes constructively with it in P waves [149]. For 
.D-waves they almost cancel each other, which seems to be consistent with the fact 
that the NN spin-orbit force is known to be very small for even orbital angular 
momentum channels. In reference [146], the uncertainties associated with spin- 
orbit terms generated by quark confinement have been emphasized by studying an 
alternative model in which confinement is described through a mass term in the 
relativistic single particle equation. The single particle Thomas term of such a model 
yields a NN spin-orbit interaction with opposite sign to the one generated by a scalar 
two-body confining potential. In summary contributions which arc usually neglected 
(antisymmetric spin-orbit terms, Thomas precession terms, etc.) might be even more 
important than those considered [65]. 

5. 5. Charge dependence of nuclear forces 

The possibility of isospin violation in the strong interaction endures as one of the 
intriguing and partially unsolved aspects in the description of nuclear forces [153]. 
After electromagnetic effects are removed from the experimental data, small but 
significant differences exist among the neutron-neutron and proton-proton interactions 
(charge symmetry breaking, CSB), and among the average of these two interactions 
and the neutron-proton one (charge independence breaking, CIB). Explicitly in terms 
of the scattering lengths, 

AacsB = a PP - a n n = 1-5 ± 0.5 fm 

1 (5-1) 
Aocib = 2 (°PP + fl nn) - a pn = 5.7 ± 0.5 fm . 

Most of the theoretical calculations to understand these numbers have been done using 
baryonic meson-exchange models that incorporate basically isospin breaking through 
the charged and neutral pion mass difference, the proton-neutron mass difference, the 
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Table 7. Different contributions to AaciB and AacsB- The number between 
parenthesis is calculated including the contribution of CSB. X indicates the 
contributions considered. 
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p° — uj electromagnetic mixing amplitude, and the electromagnetic corrections. Whilst 
the major part of the CIB has been explained by means of the mass difference between 
the neutral and charged pions, proton-neutron mass difference and p° — uj mixing have 
been used to account for most part of the CSB effects. 

From a more fundamental point of view one would expect that the source 
of Aacm and AacsB could be ultimately traced to CIB and CSB in QCD. The 
traditional quark-model approach incorporates only two sources of isospin breaking, 
namely up-down quark mass difference and the quark-exchange Coulomb interaction. 
Chiral constituent quark models incorporate two additional mechanisms contributing 
to isospin violation, the mass difference between the charged and neutral pion and 
the leading order electromagnetic loop corrections to the OGE potential (QED- 
QCD) [154]. All these contributions have been evaluated in reference [123] and the 
results are shown in tabled Taking into account the pion-mass splitting (Vqp E and 
Vqp E ), the quark-exchange Coulomb interaction (Vy) and the QED-QCD interference 
terms (Vqed-qcd) it is possible to reproduce CIB on the scattering length 
without introducing any additional parameter to those necessary for describing the 
NN interaction. The contribution of the quark-exchange Coulomb interaction and the 
QED-QCD interference terms make also possible to explain CSB on the 1 S'o scattering 
length. The influence of up-down constituent quark-mass difference is small due to 
the cancellation among different terms. 

5. 6. The NN system above the pion threshold 

Above the pion production energy the NN experimental data show several structures 
which cannot be described by simple extrapolation of models explaining the NN data 
below this energy. To correctly treat the problem it is necessary to incorporate the 
description of 7rNN states. One possible way of treating the 7rNN system is by means 
of explicit nucleon-resonance channels. The three-body problem is then reduced to a 
two-body problem assuming that the NN^ 7rNN transition only takes place trough 
resonances. Once the decay width of the resonances is incorporated in the propagators, 
the equations are able to properly couple physical two (NN) and three (7rNN) particle 
states. The effect of the coupling to resonances mainly depends on the strength 
and range of the transition potential. Therefore it is important to use a consistent 
interaction for the NN and nucleon-resonance channels. 

Experimentally it is known that the inelasticity in two-nucleon scattering up to 
at least 500 MeV in the cm. system, i.e., far beyond the two-pion threshold, is 
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(a) 



(b) 



(<-) 



Figure 23. Contributions of intermediate 7rNN states. The pion is drawn as a 
dashed line, the nucleon and A as single and double solid lines, respectively. 



predominantly due to single-pion production and it occurs essentially in the isospin- 
triplet partial waves, what suggests that pion production proceeds through single-A 
excitation. This restriction is an enormous technical simplification [155] and allows to 
model the 7rNN system by means of a NN— NA coupled channel problem. 

The main process for pion production is given by the 7rNA vertex depicted in 
figure , in which the A couples to a 7rN state. The introduction of this diagram 
has two consequences. The first one is that the OPE NA potential has a retarded 
interaction, shown in figure EJb), which modifies the AA interaction. However, the 
diagonal AA potential has a minor influence on the NN channels and usually this 
modification is not taken into account. The second one is the appearance of the delta 
self-energy diagram in the presence of a spectator nucleon shown in figure !23T c). that 
has a contribution to the A width and has to be added to the interaction. One has 
then to solve the Lippmann-Schwinger equations for the NN— NA system with the 
dressed A propagator [156], 
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Figure 25. 3 P\ NN phase shift and inelasticity. 
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E - (m A - ra N ) - - v A e (q) + ie 



(5.2) 



where q is the N— A relative momentum, fiA is the NA reduced mass and v A e (q) is 
the A self-energy interaction 

v%(q) =< 7A | -=-r-, — | 7A > (5.3) 

EA(q) - ^rr - UN + IE 

being and ldn the 7r and N energies and 7a the structure function of the 7rNA 
vertex. In equation l|5.2|) only the imaginary part of the A self-energy diagram, which 
gives the A width, is used. This is given by 



2 /^va(M m N 3 



3 Aum^ ^(feo) + m N 



(5.4) 
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Figure 26. 1 D2 NN phase shift and inelasticity. 
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Figure 27. Argand diagrams for the 1 D 2 (loft) and 3 F 3 (right) NN partial waves. 



with oj T (fco) = v rn^ + &o the pi° n energy in the 7rN center of mass system and fco the 
relative momentum of the 7rN system The effective /^(^o) can be deduced from the 
nqq vertex [157,158], 

U NA (k) = 2V2U qq F(k)e-^ (l + (5.5) 

where F(k) is the irqq vertex form factor that fixes the chiral symmetry breaking scale 
and b the wave function parameter. 

In figures [21 1251 and some NN phase shifts and inelasticities calculated in 
reference [158] are shown compared to results from the energy-dependent (dots without 
error bars in the figures) and energy-independent (dots with error bars in the figures) 
partial wave analysis of [159]. No parameters are fitted to data above the pion 
threshold. As the calculation only includes the inelasticity due to the A self-energy 
only isospin triplet partial waves are affected. For the phase shifts, the solid line shows 
results for the coupled-channel calculation including NA and AA intermediate states 
(although, as has been previously explained, the inelasticity is only included through 
the NA channel) and dashed lines refer to results including only the NN channel. For 
the inelasticities, the solid line shows the result including the energy dependent A 
width whereas dashed-dotted lines include the energy and momentum dependence, 
as discussed in reference [158]. As a general trend the experimental phase shifts are 
reasonably well reproduced by the theoretical calculation except for the P waves, 
as it happened below pion threshold. The most interesting iso-triplet partial waves 
are the 1 Z?2 and 3 F"3, because in this case the NN— NA coupling is very important. 
The Argand plot for both cases, figureEZI presents a counterclockwise behaviour with 
increasing energy, which is considered to be a signal of a possible resonance. As seen in 
figure l^Bl the resonant behaviour observed in the experimental data only appears when 
the coupling to the NA channel is considered, although the calculation overestimates 
the phase shifts as already known for models including only NA channels. 
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6. The baryon spectrum 

An exciting challenge in contemporary nuclear physics is to achieve a unified 
understanding of the baryon-baryon interaction and the baryon spectrum in terms 
of quark degrees of freedom. Once we have analyzed in detail the description of the 
baryon-baryon interaction in chiral constituent quark models with ansatz harmonic 
oscillator baryonic wave functions we should look for a consistency of such ansatz 
with the solution of the baryon spectrum from the same chiral quark-quark potential. 
Furthermore we may expect from a detailed description of the baryon spectrum to 
deepen the knowledge of some parts of the quark-quark interaction that are not very 
relevant for the baryon-baryon problem as it is the case for the confining potential. 

As we shall see these expectations will be satisfied and a quite precise description 
of the non-strange baryon spectrum will finally be achieved. In the way a long-standing 
problem in the prediction of the baryon spectrum from quark models incorporating a 
OGE plus confining quark-quark potential, say the level ordering problem involving the 
relative position of the first two excited states of the nucleon (and perhaps also of the 
A) finds a consistent solution. For the sake of clarity and simplicity and before going 
to an extensive revision of the calculational methods and results, let us establish the 
level ordering problem in the pure harmonic limit. The N*(1440) J p = 1/2+ belongs 
to the [56, + ] SU (6) Flavour- Spin x 0(3) irreducible representation and it appears in 
the N = 2 band, while the N*(1535) J p = 1/2" belongs to the [70, 1"] appearing 
in the N = 1 band. As a consequence, the N*(1440) has 2hiv energy excitation 
while the N*(1535) has only lTwo energy excitation, opposite to the order observed 
experimentally. In fact, the quadratic potential places the Roper resonance 1 GeV 
above the ground state (if its strength is fixed variationally as it is done in the RGM) , 
while it appears in the correct position (but still inverted with respect to N*(1535)) 
with a linear potential. This favors a linear form for confinement as also suggested 
by meson spectroscopy and lattice calculations [160]. We should finally say that the 
level ordering problem has been cured by means of appropriate phenomcnological 
interactions as has been the case of anharmonic terms [161], scalar three-body forces 
[162] or pseudoscalar interactions [45, 163], but without any connection to a consistent 
simultaneous description of the baryon-baryon interaction except for reference [163]. 

6.1. Formalisms of the three-quark bound state problem 

6.1.1. Feshbach-Rubinow method This method [164], applicable only for the ground 
state of the system, assumes that the ground state three-particle wave function 
depends on a single variable R = + x 2 + nx 3 ), where Xi = \ fj — r k \ with fj^ the 
particle coordinates and Cijk = 1. The variational parameter n gives account of the 
asymmetry in the force bonds and the masses. The three-body Schrodinger equation 
is then reduced to a single differential equation in the variable R with an effective 
mass and an effective potential that arc dependent on the variational parameter rj. 
The ground-state eigenvalue of the equation is minimized as a function of w. 

6.1.2. Hyper spherical harmonic method This method is specially suited for the case 
of equal quark masses. The position of three identical particles with respect to the 
center of mass can be described in terms of the Jacobi coordinates £1 and £ 2 > defined 
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in terms of the coordinates of the particles (i = 1, 2, 3) by, 
1 



6 



V2 



6 = 



(6.1) 



From them, the hyperspherical coordinates (p, 0) are defined as follows. The 
hyperradius satisfies p 2 = + £f . is a set of five angular coordinates: one 
hyperspherical angle <f> (tan0 = C1/C2) and the four angular coordinates £1 and £2 
in the three-dimensional space. 

With this choice the kinetic energy operator is 



Pi" 



-. 2 

Pi 



h 2 ( d 2 5 9 1 



TT- + TT- = hfl> + + "J^ (°) 



with 



2m 2m 
t) 2 



2m \dp 2 p dp p 2 



dcj>- 



„ + 4 cot 20 



a mi) *2(&) 



sm 



cos^ 



(6.2) 



(6.3) 



The eigenfunctions of this operator, X 2 (fi), are the orthonormalized hyperspherical 
harmonics Y[k,i 1 ,me 1 ,£ 2 ,me 2 ]- The Schrodinger equation written in hyperspherical 
coordinates becomes 

r, 2 



*(p,n) = 0(6.4) 



where V^(p, 0) is the interaction. 



6.1.3. Stochastic variational method The variational method is a widely used tool 
to solve the bound state N-particle problem. It is based on the choice of a trial wave 
function and the search of the parameter set minimizing the energy of the system. 
The trial function has to include the correlations between the different particles and 
at the same time give rise, when calculating observables, to easily calculable matrix 
elements. A linear combination of harmonic oscillator wave functions, as the simplest 

f(f i ,r j )=e- a ^-^ 2 (6.5) 

where rij refer to the particle coordinates, fulfills these two requirements and is one 
of the most popular choices in variational calculations. 

In the literature there are different strategies to select the trial wave function 
parameters. The stochastic variational method attempts to set up the most 
appropriate parameter set by following a stepwise procedure: one generates a basis 
function by choosing randomly the nonlinear gaussian parameters, ay, while the 
linear ones are determined by diagonalizing the Hamiltonian matrix. From a set of 
N randomly chosen functions one selects those which give the minimal energy. Then 
one increases the number of functions by one with randomly generated nonlinear 
parameters. Its utility is judged by the energy lost or gained and so the extra function 
is kept or discarded, respectively. This trial and error procedure is repeated until the 
basis set up leads to convergence. The versatility and the efficiency of the stochastic 
variation has been demonstrated for several systems [165]. 
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6.1. 4- Coordinate-space Faddeev method In this method the Schrodinger equation 
for the three-body problem is written as 

(Ho + V X +V 2 + V 3 ) Wi, r 2 , f%) = E^(n,f 2 ,r 3 ) (6.6) 
where Hq is the kinetic energy operator of the three particles 

H ° = t£r,- <"» 

2—1 

and Vi — Vi(fj — fife) with rjk the particles coordinates and ey^ = 1 is the potential 
energy. 

In order to guarantee a unique solution of the three-body problem Faddeev 
decomposed the wave function of the system as [166] 

* = fa + fa + fa , (6.8) 

where the Faddeev components satisfy the coupled equations 

(E-Ho-Vi)<p i = V i (ct> j + <Pk) ■ (6.9) 

This is the form used in coordinate-space calculations where the spatial part of fa is 
decomposed in terms of the spherical harmonics relative to the Jacobi coordinates as 
defined in equation l|6.1fl . 

Due to the presence of confining and constant terms in the potential, the 
convergence of these equations with respect to the number of angular momentum 
channels is rather slow. Therefore, a modification of the method has been proposed in 
references [167,168] in order to accelerate the convergence of the solution with respect 
to the number of angular momentum channels. The two-body interaction is splittcd 
into confining and non-confining parts as, 

Vi = V? + Vf c , (6.10) 

which leads to the modified form of the Faddeev equations 

(E-H Q - Vf - V 2 C - V 3 C - V t NC ) fa = V» c (fa + fa) . (6.11) 

This modified method, however, must be used with some care since the splitting 
into confining and non-confining parts is not unique so that it can lead to spurious 
solutions [167,168]. 

6.1.5. Momentum- space Faddeev method Formally equation Ij6.9|l can be rewritten 
as 

fa = G U(fa+fa) (6.12) 
where Go is the propagator of three free particles 

Go^(E-Ho)- 1 (6.13) 
and ti is the two-body T-matrix which obeys the Lippmann-Schwinger equation, 

U =V + V t G U (6.14) 

The method to solve the Faddeev equations for three quarks in momentum space 
has been described in detail in reference [169]. Let us mention here that there is a 
problem when dealing with confining potentials of the form 

V c (r)=br n ; n = l,2,... (6.15) 
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in momentum space, since the Fourier transform does not exist because Vc — > oo 
when r — > oo. However, if the potential I|tj.l5fl is replaced by the finite potential 

t/y \ (b(r n -R n ) ;r<R . . 

V(r) = { Q , r>R (6.16) 

for which the Fourier transform is well defined, in the limit R — > oo both potentials 
will give the same result provided that the eigenvalues are related by 

E c — E — 3bR n (6.17) 

where Ec and E are the eigenvalues corresponding to the interactions Ij6.15|l and 
H6.16JI . respectively. 

If there are no tensor or spin-orbit forces the Faddeev equations for the bound- 
state problem of three quarks can be written as 

1 1 f 1 f°° 

< Pm ; £iXiSiTi\<pf ST >= _ E 2 ^ L q > d<!j 

x tt s < T <( Pi ,p>;E- q y2 Vi )Af^{ P ' iq%Pjqj ) 

x < SiTtlSjTj > ST < P3 q 3 ;l^ J S ] T ] \4> I ; ST >, (6.18) 

where Si and Ti are the spin and isospin of the pair jk while S and T are the total 
spin and isospin. £i ( P i) is the orbital angular momentum (momentum) of the pair 
jk, \i (qi) is the orbital angular momentum (momentum) of particle i with respect to 
the pair jk, and L is the total orbital angular momentum. cos8 — ql ■ qj/{qiqj) while 

m j m k , Rin x 

m = — , (6-19) 

m,j + m k 

= mi {m 1+ m k ) ^ 
rrii + rrij + m k 

are the usual reduced masses. For a given set of values of LST the integral 
equations l|6.18|l couple the amplitudes of the different configurations {^AiS 1 ^} with 
(— yi+ s t+ T i — 1 as required by the Pauli principle since the wave function is colour 
antisymmetric. 

The spin-isospin recoupling coefficients < SiTi\SjTj >st are given by 



< SiTilSjTj > S T= (-) 5 *+°*- s \fi2Si + l)(25j + ijwiajVkSaaSiSj) 

x{-) T ^- T ^/{2T l + 1)(2T, + 1) W(TjTkTTf, TiTj), (6.21) 

with a \ and r, the spin and isospin of particle i, and W is the Racah coefficient. The 
orbital angular momentum recoupling coefficients A£ * J 3 ' ( P iQi P jQj) are given by 

itli Xilj Xj / i \ ' (-iliXiL ^itjXjL V V V 

A L xPiliPjlj) ~ 2L + X mi,M—mi,M mj,M—mj,M XtM—mi* Ijirij 

xrx jM - mj cos[-M(qj,ql) - m l {q i , P ' i ) + mj(qj, Pj )], (6.22) 
with T£ rn = if £ — m is odd and 



(-)(<+")/y(2l+l)(l + m)!(l-m-)I 

2 e ((£ + m)/2y.((£-m)/2)] [ ' ' 
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if £ — m is even. The angles (qj,<li), (<fi,.p/), and (qj,Pj) can be obtained in terms of 
the magnitudes of the momenta by using the relations 

P'i = ~Qj - (6.24) 

&=?i+^ L & (6-25) 
m k 

where ij is a cyclic pair. The magnitude of the momenta p\ and pj , on the other hand, 
are obtained in terms of qi, qj, and cosO using equations 1)6.24(1 and ((6.25(1 as 



P'i = \M+ (—Xlf + —Mjcose, (6.26) 
V J \m k l m k 



,2 I '// \ ,,J i 



Pj = \ii+ — qj + —q l q 3 cos6. (6.27) 

y \m k J J m k 

Finally, the two-body amplitudes t i Ti {pi,p'f, E — qf/2i>i) are given by the 
solution of the Lippmann-Schwinger equation 

/>oo 

#**ME-tf/toi)=V***M+ pfdp'lV^^p'!) 

Jo 

with 

V- iSiTi M) = - rldriuMnWf^ir^uMn). (6.29) 

JO 

and j£ the spherical Bessel function. 



6. 2. Results of different potential models 

6.2.1. Isgur-Karl model In a series of pioneering works Isgur and Karl [19] performed 
an extensive study of the baryon spectrum. They used an interaction composed of 
a confining potential, usually taken to be a harmonic oscillator, the hyperfine term 
of the OGE potential proposed by de Rujula, Georgi and Glashow [16], and finally 
anharmonic terms to treat on a different footing states belonging to different major 
shells. They worked in first-order perturbation theory, diagonalizing the Hamiltonian 
in a harmonic oscillator basis, and they obtained an overall good agreement with 
experimental data. Nonetheless to get such agreement the zero order masses of 
positive- and negative-parity states were fitted independently, without any sound 
justification. 



6.2.2. Bhaduri-like potential models Bhaduri et al [20] performed a calculation 
of the ground-state baryon masses using an interaction designed to describe the 
meson spectra, composed of a standard OGE interaction besides a linear confinement. 
They used the Feshbach-Rubinow method assuming only 5-wave interactions. They 
obtained M N = 1052 MeV and Ma = 1354 MeV, which are more than 100 MeV 
above the experimental values. 

Later on a complete coordinate-space Faddeev method calculation was carried 
out with the same interaction [21]. The comparison of the results with those of the 
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harmonic oscillator diagonalization method allowed to establish the high degree of 
accuracy of the last technique. For the masses of the N and A the values 1024 MeV 
and 1318 MeV, respectively, were predicted. These are about 80 MeV above the 
experimental values although their splitting is correctly given around 290 MeV. For 
the N and A excited spectra (see figure |2EJl they got a reasonable description of the 
negative parity excitations while in the positive parity case the first radial excitation 
for both, N and A, was above the data by about 300 MeV. 

The problem of the relative position of the negative- and positive-parity 
excitations of the nucleon was addressed in reference [162] by implementing the 
Bhaduri potential with a scalar three-body force corresponding to the exchange of 
two cr's. Hyperspherical harmonic and coordinate-space Faddeev calculations were 
consistently performed. This allowed for a resolution of the level ordering problem in 
both cases, N and A, and for a pretty good description of the non-strange baryon 
spectra. The question immediately arising was whether this effective three-body 
force effect could be obtained or not with a two-body description once the chiral 
pseudoscalar exchanges were incorporated in the quark-quark potential. 



6.2.3. Goldstone-boson exchange potential model In this model the interaction 
between quarks is assumed to be given by a linear potential responsible for confinement 
plus a one-boson-exchange potential generated by the exchange of the octet of low- 
mass pseudoscalar mesons and a flavour singlet 77' [45,170], i.e., 

V = V x + V conf (6.30) 



where 



v° ctet (nj) 



EMr«)A?AJ + 7,(r5;-)AfA 



a=4 



V° in9let (r7 j ) = -V v ,<j i 



(6.31) 
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Figure 29. IJeldl^ — \ SjT-and sp^lriga. of reference [170]. 



The model was first applied by Glozman and Riska [45] within a first-order 
perturbation theory approach. In the more elaborated calculations of the model 
the one-boson-exchange potentials were constructed with different meson-quark-quark 
coupling constants for the octet and singlet and were regularized at the origin by 
smearing the ^-function part with appropriate cut-off parameters which are different 
for each meson. As the coordinate-space Faddeev method allows to obtain a solution 
only in the case of a non-relativistic kinetic energy operator and the effect of replacing 
the non-relativistic kinetic energy operator l|6.7|l by the relativistic one, i.e., 



(6.32) 



was estimated to be important for the level ordering problem, a calculation in the 
stochastic variational formalism was carried out (see figure |2*§1 from [170]). The fit 
to the experimental spectrum was remarkable up to 900 MeV excitation energy. The 
first positive- and negative-parity nucleon states were correctly predicted though some 
caution should be in order to judge the A case. Actually the Roper of the A was 
predicted above the first negative-parity state though experimental error bars do not 
allow in this case to draw any definitive conclusion. 



6.2.4- Chiral constituent quark models The non-strange baryon spectrum has been 
also studied from the chiral quark potential of section Reference [44] presented a 
calculation with a small value for ro (ro = 0.0367 fm), the regularization parameter 
for the contact term of the OGE, and truncating the Hilbert space by means of 
a hypersphcrical harmonic approach up to k = 2 excitations. Similar results were 
presented in reference [171]. These results were strongly criticized [172] by arguing 
that such a small value of ro would give rise to unstable values for the baryon masses. 
Moreover, it was argued that the level ordering problem could not be solved with 
a chiral constituent quark model due to the presence of the OGE term and the 
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Figure 30. N(l/2+) and A(3/2+) ground state masses as a function of ro. 



subsequent reduction of the chiral quark coupling with respect to the Goldstone-boson 
exchange potential case. However the results of [44] and [171] can be alternatively 
reproduced almost exactly from a momentum-space Faddeev calculation including 
only the lowest order configurations (£i, Aj, Si,Ti) [163,173] for a much bigger value 
of ro = 0.2 fm, lying in the stability region. The dependence of the N(l/2 + ) and 
A(3/2 + ) ground state masses on the regularization parameter ro, as given by the 
Faddeev calculation, is illustrated in figure [3U1 While the baryon-baryon interaction 
does not crucially depend on ro, the baryon spectrum does. For values of ro below 
0.1 fm the N(l/2 + ) mass decreases very quickly while the A(3/2 + ) remains almost 
stable. As a consequence, the excitation energy of the whole spectrum would become 
too high. This result can be easily understood taken into account that the short- 
ranged spin-spin force is attractive for the nucleoli but repulsive for the A, lowering 
the L — nucleon states. It is also worth to notice that the nucleon ground state 
energy almost does not change for ro values greater than 0.1 fm. 

Increasing the value of ro the excitation energy of all the states diminishes, but 
maintaining a similar ordering. The excitation energy can be easily increased just 
by changing the confinement strength and slightly modifying the strong coupling 
constant. The spectrum obtained in a lowest-order Faddeev calculation with the set of 
parameters A of table is shown in figure |^ where the solid lines correspond to the 
predicted excitation energies with respect to the nucleon ground state and the shaded 



Table 8. Different sets of quark model parameters. 



Set 


A 


B 


C 


D 


m 9 (MeV) 


313 


313 


313 


313 




0.72 


0.65 


0.485 


0.50 


a c (MeV ■ fni" 1 ) 


72.518 


60.12 


67.0 


110.0 


OLch 


0.0269 


0.0269 


0.0269 


0.0269 


r (fm) 


0.2 


0.8 


0.25 


0.74 


m a {im~ 1 ) 


3.42 


3.42 


3.42 


3.42 


m 7r (fm — 1 ) 


0.7 


0.7 


0.7 


0.7 


A^fm" 1 ) 


4.2 


5.4 


4.2 


2.0 


A CT (fm- 1 ) 


4.2 


4.2 


4.2 


2.0 
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Figure 31. Relative energy N and A spectrum up to 1 GeV excitation energy 
for the set of parameters A (left) and B (right) of table 151 



Table 9. Convergence of baryon masses (in MeV) with respect to the number 
(Nr.) of configurations (£}, \j, Si,Ti). See text for details. 



Nr. 


M N 


Nr. 


M A 


Nr. 


-^N*(1440) 


Nr. 


A^N*(1535) 


2 


115 


1 


423 


2 


641 


4 


643 


■1 


68 


2 


419 


4 


622 


8 


601 


6 


21 


3 


418 


6 


587 


12 


570 


8 


12 






8 


578 


16 


551 


10 


2 






10 


562 


20 


540 


12 









12 


558 







areas to the experimental values including their uncertainties. As can be seen, apart 
from the energy difference between the positive and negative parity excitations of the 
nucleon, the spectrum is almost the same as those reported in [44] and [171]. These 
results make clear that the regularization depends on the model space in which the 
calculation is done and the parameter of this regularization should not be understood 
as a true parameter of the model Hamiltonian but instead as a free parameter to be 
fitted. 

Let us note that the level ordering problem remains in spite of the fact that the 
(c?i • dj ) (fi ■ Tj ) structure of the pseudoscalar potential tends to favor the right location. 
Indeed, the correct ordering of the Roper and the negative parity excitations of the 
nucleon can be obtained by increasing the strength of the pseudoscalar potential. The 
strong correlation among the parameters of the model obliges for a modification of 
a s and a c obtaining the set of parameters B of table |H1 The spectrum obtained with 
this set of parameters is shown by solid lines in figure El It can be checked how 
the correct level ordering between the positive and negative parity can be obtained at 
the expense of having a stronger OPE (the dashed line in this figure stands for the 
results of reference [44]). However, enhancing artificially the pseudoscalar potential 
the agreement for the two-baryon sector, specially the binding energy of the deuteron, 
is destroyed [46] . Therefore in the lowest-order Faddeev calculation the level ordering 
problem has no solution if consistency with the two-baryon system is required. 

Let us now examine the situation when one extends the calculation to include 
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all the configurations (^j, Aj, Si, Tj) with and up to 5. The convergence with 
respect to the number of configurations is shown in table ED for N, A, N*(1440), and 
N*(1535) states corresponding to the set of parameters B of table IS1 (figure l?T|) . The 
mass difference with respect to the nucleon ground state mass as the reference (i.e., 
Mn = 0) is given. As can be seen, the convergence with respect to the number of 
configurations is different for different states. So, the A comes down by only 5 MeV 
while the N, N*(1440), and N*(1535) come down by approximately 100 MeV when 
one includes the higher orbital angular momentum configurations. It is also worth to 
notice that while in the lowest-order calculation the N*(1440) lies below the N*(1535), 
when one includes the higher orbital angular momenta the situation is reversed. 

Thus, the higher orbital angular momentum components play an important role 
on the description of the low-energy baryon spectrum. They influence the relative 
energy of the nucleonic and A sectors and they also affect the relative position of 
the positive- and negative-parity excitations of the nucleon. In figure [221 we show 
the baryon spectrum obtained with an extended Faddeev calculation for the set of 
parameters C of table |S1 that correspond to the ones used for the description of the 
baryon-baryon interaction [43,62]. The parameter rg has been taken as 0.25 fm, 
a typical value for spectroscopic models, such that the delta function of the OGE 
interaction is regularized in the stability region (see figure 001 • As one can see, the 
relative position of the positive-parity states N, A, and N*(1440) are correctly given, 
although the negative-parity states lie below the experimental results. It is interesting 
to check that the spectrum obtained is qualitatively very similar to the one obtained 
in [44]. The slope of the confining potential, a c , which in set C of tablets 67.0 MeV 
fm -1 , is not so different from the value 72.5 MeV fm -1 used in [44] ft. 

It is worthwhile to investigate the role played by the OGE and OPE potentials 
in the level ordering problem. In reference [163] the energies of the N*(1440) and 
the N*(1535) were calculated starting with the set of parameters C and varying the 
contribution of the colour-magnetic interaction by modifying the coupling constant a s 

(I One should not forget at this point that consistency with the NN sector does not impose any 
restriction to the value of a c , because the confining interaction does not contribute appreciably to 
the NN potential 
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Figure 32. Relative energy N and A spectra up to 1 GeV excitation energy 
including all configurations with I and A up to 5 for the set of parameters C. 
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Figure 33. N* (1440) and N* (1535) masses as a function of the OPE cutoff mass. 

of the OGE interaction. It turns out that the relative position of the N*(1440) and the 
N*(1535) is not very much affected by the modification of a s - This behaviour indicates 
that the correct level ordering between the negative and positive parity excited states 
of the nucleon does not come from a weakening of the OGE interaction. Furthermore, 
a suppression of the OGE would imply a stronger pseudoscalar interaction in order 
to reproduce the A— N mass difference, and therefore one would obtain a model that 
it is incompatible with the understanding of the basic features of the two-nucleon 
system [46]. 

Regarding the pseudoscalar interaction the (a ■ a) (r • f) structure of the OPE 
gives attraction for symmetric spin-isospin pairs and repulsion for antisymmetric ones 
(a quite distinctive feature since the colour-magnetic part of the OGE gives similar 
contributions in both cases). This lowers the position of the first nucleon radial 
excitation, N*(1440), with regard to the first negative parity state solving part of 
the discrepancy between usual two-body potential models. To illustrate this point 
the energy of the N*(1440) and the N*(1535) was recalculated starting with the set of 
parameters C of tabled but increasing the contribution of the pseudoscalar interaction 
by letting the OPE cutoff parameter A v to increase. The results are shown in figureOSl 
As can be seen, the inversion of the ordering between the positive and negative parity 
states can be achieved if A w becomes sufficiently large (around 7 fin -1 ). A model with 
such a strong cutoff for the OPE is not realistic, the resulting A— N mass difference 
being around 955 MeV. If one fits again the A— N mass difference by modifying the 
confinement constant and suppressing the OGE one loses again the inversion between 
the negative and positive parity states. Therefore we are forced to conclude that a 
consistent solution to the level ordering problem cannot be attained with the pure 
non-relativistic chiral constituent quark model. 

Recently the baryon spectrum has been studied by means of a chiral constituent 
quark model incorporating relativistic kinematics [49] . We show in figure 02] the 
spectrum obtained from a momentum-space Faddeev calculation, in order to include 
full three-body relativistic kinematics using the set of parameters given in column D 
of table [8] As one sees the level-ordering problem that was evident in the spectrum of 
figure |23 has now been cured. As suggested long-time ago the relativistic kinematics 
has an important influence on the N*(1440) mass [174, 175] and the correct level 
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Figure 34. Relative energy N and A spectra up to 1 GeV excitation energy with 
relativistic kinematics for the set of parameters D of table IS1 

ordering is obtained due to the combined effect of relativistic kinematics and one- 
pion-exchange. The OGE does not play a major role to this respect although it 
should be remarked that a decreasing of the strength of the OGE goes against the 
correct level ordering. 

From a technical point of view, the replacement of the non-relativistic by the 
relativistic kinematics for short-range attractive interactions has to be done together 
with a careful study of the short-range part of the potential in the relativistic approach. 
This is due to the fact that the non-relativistic propagator falls down as 1/k 2 while the 
relativistic one does as 1/k. This is, for example, the reason why the semirelativistic 
solutions of reference [176] collapse when a scalar and an one-gluon-exchange potential 
are added to the pseudoscalar interaction. This makes evident the risk of replacing 
the non-relativistic kinetic energy by the relativistic one when the details of the short- 
range part of the potential are not well known. 

6.2.5. Other potential models In view of the importance of relativity with regard 
to the inversion of positive- and negative-parity nucleon states it is also interesting 
to mention studies using a fully covariant approach [177-180]. This approach is 
based on the assumption of instantaneous interactions between the quarks which has 
the consequence that the three-quark Bethe-Salpeter equation reduces from its four- 
dimensional form to a three-dimensional form, the Salpeter equation, whose degree of 
complexity is similar to the non-relativistic Faddeev method. 

A very complete study of the baryon spectrum has been carried out by the Bonn 
group [178-180] within a covariant formalism using the Salpeter equation and where 
the spin has been treated dynamically. Its model consists of a constant term, a linear 
confining term, and a residual quark-quark interaction. Two models of the part of 
the interaction containing the constant and confining terms (models A and B) which 
correspond to different choices for the Dirac structure of each of these terms are 
considered. The residual quark-quark interaction is based in 't Hooft's instanton- 
induced interaction. This determines its Dirac structure while its radial structure is 

I "* 2 / X 2 

taken of gaussian type e -|rij ' with A = 0.4 fm. This instanton-induced interaction 
plays the role of the one-gluon-exchange potential in the more standard models. Since 
this interaction acts only in flavour-antisymmetric states it is responsible for the 
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Figure 35. Main hypcrradial component of the nucleon wave function and 
diagonal kernels of the NN potential V(R) for the (S, T) = (1,0) channel. The 
solid line corresponds to the solution of the Schodinger equation, the long, medium 
and short-dashed lines to gaussian ansatz with values of the parameter b = 0.4, 
0.518, and 0.6 fm, respectively. 

hyperfine splitting in the baryon spectrum. Thus, the A spectrum where all the states 
are flavour-symmetric is determined only by the constant and linear confining terms. 
This fact is used to fit the parameters of the two models (A and B). The strength of 
the instanton-induced residual interaction is then fitted to the A— N mass splitting. 
A reasonable description of the A spectrum is obtained for both models. However, in 
the case of the nucleon spectrum only model A gives a satisfactory description, the 
problem of the inversion of positive- and negative-parity nucleon excitations being only 
partially cured since complete inversion is not achieved. This shows that relativity 
alone is not sufficient to give rise to the inversion of states observed in the experimental 
spectrum (as also demonstrated in reference [169]) being mandatory to this aim the 
inclusion of a pseudoscalar exchange interaction [49] . 

Let us finally mention that in reference [181] the mass spectrum of the low-lying 
ground-state baryons has been calculated by solving the Dirac equation for each quark 
in a single-particle confining potential and treating the residual interactions (OGE and 
OPE potentials) as low-order perturbations. 

6.2.6. Consistency of the baryon spectrum wave functions with the ansatz baryon- 
baryon wave functions 

After discussing the one-baryon problem one question should immediately arise. 
The two-baryon system has been studied by means of gaussian-like wave functions for 
the spatial part while later on the one-baryon problem has been solved and therefore 
one has the baryon wave functions predicted by the quark-quark Hamiltonian. Does 
it make sense the gaussian hypothesis for the radial part of the wave function? To 
answer this question figure 1851 compares the wave function predicted by the quark- 
quark Hamiltonian used for figure 1^21 to the ansatz used to study the baryon-baryon 
interaction for three different values of the oscillator parameter. Although at long 
range all the wave functions look very similar, the discrepancy at the origin may reach 
a factor four. This effect may be very important for quantities depending on the value 
of the wave function at small distances. Hence an evaluation of the importance of this 
difference in the calculation of the phase shifts should be done. From the technical 
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point of view a complete RGM calculation would be very difficult to perform with 
the full solution of the Schrodinger equation, but one can easily compare the diagonal 
kernels of the NN potential (which correspond to the diagonal part used to calculate 
the phase shifts and the deuteron properties) with both type of wave functions as 
shown in figure |23 We see how the NN potential calculated with the full solution 
of the Schrodinger equation looks pretty similar to the potential calculated with the 
harmonic oscillator wave function with b =0.518 fm, this is precisely the value used 
in the RGM calculation of the phase shifts and deuteron properties conferring to it a 
self-consistent character. Other values of b give quite different potentials even in the 
medium- long range part. In this way the validity of the gaussian ansatz wave function 
in the RGM calculations is confirmed, and the old controversy about the possible 
values of the parameter b can be solved. 

6.3. Exotic baryons 

Once the baryon spectrum has been analyzed one may wonder about the predictions 
of chiral quark models for exotic baryons, i.e., systems with baryon number +1 whose 
quantum numbers cannot correspond to a simple three-quark picture. There has been 
a renewed interest on this subject from the announcement of the discovery of the so- 
called + resonance, with strangeness +1, in several experiments listed in Table ITUl 
This resonance, carrying unit positive charge, shows out as a peak in the invariant 
mass distribution of (K + n) or (Kgp). The location of the peak corresponds to an 
invariant mass that ranges from 1520 to 1570 MeV depending on the experiment and 
a small width in all cases lesser than 25 MeV (actually KN phase shifts analyses 
give an upper bound of ~ 1 MeV, see reference [208] and references therein). The 
angular momentum and parity of the resonance are still experimentally undetermined 
whereas there are experimental indications on its iso-singlet character [186]. After 
the initial excitement several recent experiments also listed in Table HUI have failed to 
find evidence for pentaquark signals (see Ref. [209] for a detailed description of the 
experiments). 

Different possible explanations for such a peak have been discussed. Leaving aside 
its possible origin as a kinematic reflection of the intermediate production of known 
meson resonances [210] what deserves a complete analysis and having been discarded 
its KN scattering potential nature, the most compelling explanations are associated 
to a chiral soliton or a pentaquark (ududs) structure. The chiral soliton model by 
Diakonov, Petrov and Polyakov [10] predicts an exotic l/2 + baryon belonging to 
the 10 f representation of SU (3) Flavour with a mass around 1540 MeV and a width 
less than 15 MeV (this number has been recalculated as 30 MeV in reference [211]) 
stimulated its experimental search. It should be also mentioned that years before 
Praszalowicz [212] predicted with a Skyrme model a 1530 MeV mass for the iso-singlet 
member of the antidecuplet. 

From a quark model point of view and under the assumption of a J p = l/2 + 
assignment, alternative mechanisms, mainly based on colour or chiral interactions, 
have been proposed. 

Among the colour based models we shall mention the diquark-diquark and 
diquark-triquark models. One should first realize that a simple independent quark 
model based on a confinement plus OGE potential would predict the ground 
state pentaquark to have negative parity, since the four quarks and the antiquark 
would be in 0s states and the intrinsic parity of the antiquark is negative. This 
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Table 10. Overview of the positive and null pentaquark experiments. 

Positive experiments 



Exp. 



Reaction 



Mass (MeV/c 2 ) Width (MeV) Significance Rcf. 



LEPS(l) 

LEPS(2) 

DIANA 

CLAS-d 

SAPHIR 

CLAS-p 

ITEP 

SVD 

HERMES 

ZEUS 

COSY 



7 12 C -> K+K-X 
~/d^ K+K-X 
K+Xe -» X s °p(Xe)' 
jd - 
IP - 
IP - 
vA- 
pA - 
e+d 
ep — 
PP - 



■ K+K-p(n) 

■ K+K° s (n) 

■ K+K-n+(n) 

* K° sP X 

♦ P K° S X 
-» K° s pX 

K° sP X 



1540±10 

- 1540 

1539±2 

1542±5 

1540±6 

1555±10 

1533±5 

1526±3 

1526±3 

1522±3 

1530±5 



< 25 

< 9 

< 21 

< 25 

< 26 

< 20 

< 24 
13±9 
8±4 

< 18 



4.6a 

4.4(7 

5.2cr 
4.8cr 
7.8a 
6.7a 
5.6cr 

<~ 5(7 
<~ 5(7 
~ 5(7 



[182] 
[183] 
[184] 
[185] 
[186] 
[187] 
[188] 
[189] 
[190] 
[191] 
[192] 



Null experiments 



Exp. 



Reaction 



Rcf. 



HERA-B 

BES 

CDF 

BaBar 

ALEPH 

DELPHI 

FNAL E690 

HyperCP 

PHENIX 

BELLE 

FOCUS 

LASS 

L3 

SELEX 
SPHINX 



pA -» e+x 

e+e- -» J/V#(2S)] -» ee 

pp -» e+x 

e+e- -» T(45) 
e+e~ — > Z — > 
e+e - —fZ—tqq 

PP ~ > pK~TT+Q+ 

(tt+,K+, P )Cu-> G+X 
dAu -» O+A 

ATA -» e+x 
7 p -» e+x 

A+p -» X+mr+ 

77^ ee 

(7r,p,E)p^e+X 
pC(JV) -» e+A°C(A) 



[193] 
[194] 
[195] 
[196] 
[197] 
[198] 
[199] 
[200] 
[201] 
[202] 
[203] 
[204] 
[205] 
[206] 
[207] 



difficulty can be overcome by introducing multiquark correlations. Jaffe and 
Wilczck [213], realizing the presence of strong colour-spin correlations in the baryon 
spectrum, built a diquark-diquark-antiquark ((ud)(ud)s) model with the diquark 
structures {^Flavour, ^Colour, S — 0) favoured by the colour magnetic (spin-spin) OGE 
interaction. The assignment of the nucleon Roper resonance as well as the to the 
same 10f representation allowed for an understanding not only of the pentaquark 
mass but also of the Roper mass on the basis of three-quark 8_p— pentaquark 10f 
mixing. On the other hand Karliner and Lipkin [214] proposed a diquark-triquark 
model ((ud)(uds)) with a triquark (6f,3c,S = 1/2) structure so that the diquark 
and triquark have a colour electric OGE interaction (the Coulomb term) between 
them which is identical to the quark-antiquark interaction in a meson. Both models 
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give the right mass for the pentaquark. The small width is attributed to the weak 
coupling to the KN continuum due to the difference in colour, spin and spatial wave 
functions. Both models predict new stable pentaquarks in the heavy quark sectors 
what could serve, if experimentally found, to disentangle the role of quark correlations. 

On the other hand chiral quark models may provide a natural explanation for the 
positive parity ground state pentaquark. This is not surprising since they incorporate 
Goldstone bosons merging from spontaneous chiral symmetry breaking and it is 
precisely this breaking the crucial fact in the formulation of the chiral soliton model. 
Explicitly the flavour-spin interaction may shift the orbitally excited four-quark state 
below the (Os) 4 since the L = 1 configuration allows for a completely symmetric 
flavour-spin four quark wave function whereas the (Os) 4 configuration does not [215]. 
This parity inversion is similar to the one giving rise to the correct level ordering in 
the nucleon spectrum as we have seen. Stancu and Riska [216] have pointed out that 
pentaquark stability could be provided by a strong spin-spin interaction between the 
light flavour and the strange antiquark. 

These arguments should encourage more work along these lines in order to get a 
deeper understanding at the theoretical level. From the experimental point of view 
further experiments are needed to clarify the situation. 
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7. Dibaryons and tribaryons 

The possible existence of unstable non-strange two- and three-baryon resonances 
corresponding to bound-state solutions of systems composed of nucleons and deltas 
has been speculated for years [217-226]. In the case of two-body systems (that we 
shall refer to as dibaryons) they will decay mainly into two nucleons and either one or 
two pions, while for the three-body case (tribaryons) they will decay mainly into three 
nucleons and either one, two, or three pions. In principle, any nucleus with at least 
three nucleons can serve as test system that may be excited by forming a dibaryon or 
a tribaryon [227]. The lifetime of bound states involving one or more A's should be 
similar to that of the A in the case of very weakly bound systems and larger if the 
system is very strongly bound. Therefore, these resonances will have widths similar 
or smaller than the width of the A so that, in principle, they may be experimentally 
observable. 

7.1. Dibaryons in the constituent quark model 

The study of dibaryon resonances within the framework of the constituent quark model 
has generally been carried out by considering the baryon-baryon interaction obtained 
from a model consisting of confinement and a OGE interaction. This baryon-baryon 
interaction exhibits a repulsive core at short distances in most channels due to quark 
antisymmetrization effects and the hyperfine quark-quark interaction. Thus, dibaryons 
may exist in such channels where the strong short-range repulsion is missing and 
therefore quark-exchange forces can produce a tight deep bound state of six quarks. 

The detailed investigation of the nature (attractive or repulsive) of the different 
channels of the NN, NA, and AA systems has been an important subject of study 
ever since the first applications of the constituent quark model to derive a baryon- 
baryon interaction. Oka and Yazaki [30] showed that the Pauli principle and the 
OGE interaction generate a hard core in the NN (j, i) (j stands for the total angular 
momentum and i for the isospin) channels (0,1) and (1,0), whereas in the AA channels 
(0,3) and (3,0) the interaction is strongly attractive such that the (j, i) = (3, 0) channel 
is bound. In a later revised analysis by Oka [95] the AA (j, i) channels (1,0), (0,1), 
and (3,0) were found to be attractive while the channel (0,3) was weakly repulsive. 

The AA states = (0,3) and (3,0) have also been studied by Maltman [218] 
using the standard confinement and OGE interaction model constructed in [228]. He 
found that the (3,0) state is bound by 260 MeV while the (0,3) state is bound by only 
30 MeV which is in qualitative agreement with the results of Oka and Yazaki [30,95] 
and Cvetic et al [98]. 

The AA state (j, i) = (3, 0) has been called the deltaron by Goldman and 
collaborators [96,97] since this state is the analogue of the deuteron in the case of two 
A's. Their model incorporates the OGE interaction as well as the quark derealization 
(similar to electron derealization in nuclear physics) and color screening (predicted 
by unquenched lattice calculations). Indeed the combined effect of the last two terms 
produces the intermediate-range attraction observed in the NN system. The binding 
energy of the (3,0) state in their model is larger than 200 MeV so that the state lies 
below the 7rNA threshold. 
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7.2. Dibaryons and tribaryons in the chiral constituent quark model 

This study of the possible existence of such structures has been done with the help 
of the baryon-baryon interactions obtained from the chiral constituent quark model 
described in section [5J It presents the advantage that once the model parameters are 
fixed in the study of the NN system, all the other interacting potentials are parameter- 
free. Besides, it allows for a comparison of the results obtained by means of the 
local (BO potential) and non-local (RGM potential) interactions derived from the 
underlying quark-quark potential. 

In order to perform the calculations one can assume that the A is a stable particle, 
that is, one can neglect its width and the effects of retardation in the OPE interaction 
of the baryon-A subsystems. These two effects have been estimated in the case of the 
simple NA system [229] . There, it was found that the assumption of a stable A leads to 
very reliable predictions for the mass of NA resonances since the effects of retardation 
and width of the A are responsible for producing the width of the NA resonance but 
have almost no effect over its mass. Thus, one can be confident that the masses of 
the bound states containing A's will not change very much when its unstable nature 
is explicitly taken into account. For the three-baryon systems one can take advantage 
of the experience gained in the three-nucleon bound-state problem [230,231], where 
one knows that the dominant configuration is that in which all particles are in S-wave 
states. However, in order to get reasonable results for the binding energy, the S-wave 
two-body amplitudes used as input in the Faddeev equations must also contain the 
effect of the tensor force. So, for example, in the case of the Reid soft-core potential if 
one considers only the S-wave configurations but neglects the tensor force in the two- 
body subsystems there is no bound triton. However, including the effect of the tensor 
force in the NN 3 Si — 3 Di channel, but using only the x Sq and 3 Si components of the 
two-body amplitudes in the three-body equations (2-channel calculation), one gets a 
triton binding energy of 6.58 MeV against the experimental value 8.49 MeV. Notice 
that including the remaining configurations (34-channel calculation) a triton binding 
energy of 7.35 MeV is obtained [232]. This means that the S'-wave truncated T-matrix 
approximation leads to a binding energy which differs from the exact value by less than 
1 MeV. Hence this approach may be efficient to look for the best candidates for bound 
states and the energy ordering of the different tribaryon systems. 



7.2.1. Formalism of two-body bound states Let us consider two baryons Bi and Bi 
in a relative S'-wave interacting through a potential that contains a tensor force. Then 
there is a coupling to the B1B2 D-wave so that the Lippmann-Schwinger equation of 
the system is of the form 

I'.s'/ 



E-p[ 2 /2 m +it ^ 

where ji and %i are the total angular momentum and isospin of the system, while ijSj, 
X'iS 1 ^ and l"s" are the initial, intermediate, and final orbital angular momentum and 
spin of the system, respectively, pi and rji are, respectively, the relative momentum 
and reduced mass of the two-body system. Tables El and El show the corresponding 
NN, NA and AA two-body channels in a relative S-wave that are coupled together for 
the possible values of j and i (since the NN state is the one with the lowest mass, in 
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Table 11. NN channels (^nn>Snn)i NA channels (1nAi s na)i an d AA channels 
(ZaAi^Aa) that are coupled together in the 3 Si — 3 Di and 1 So NN states. 



NN state j i (£nn,snn) 


(Ina, «na) 


(^AA , s Aa) 


3 5i- 3 Di 1 (0,1), (2,1) 




(0,1),(2,1),(2,3) 


1 S'o 1 (0,0) 


(2,2) 





this case we have considered also the possibility of transitions to larger mass states like 
NA and AA). In the NN and AA cases the Pauli principle requires (— y<+ 5 *+'< = — 1. 

As mentioned before, for the solution of the three-body system one can use only 
the component of the T-matrix obtained from the equation (|7.1(l with U = l" = 0, so 
that for that purpose one can define the 5*-wave truncated amplitude 

t i , SiU {p h p'!-,E)^t^;{jp u phE). (7.2) 

7.2.2. Formalism of three-body bound states If one restricts to configurations where 
all three particles are in S'-wave states, the Faddeev equations for the bound-state 
problem in the case of three particles with total spin S and isospin I are (see equation 

i r°° c 1 

T i$i(j>i1i) = y^y^< s Al s j«j >sj 2 / 9jdqj / dcosdti. Siii (pup^; E - ^ 2 /2^) 

where pi and qi are the usual momentum-space Jacobi coordinates (see section IB*. 1.5(1 . 
r\i and i/, are the reduced masses defined by equations (|6.19() and l|6.20|l . and 
< Siii\sjij >si are the spin-isospin coefficients given by equation i|6.21|l . 

The three amplitudes T*£}(piqi), T^ 2 I (p 2 q 2 ), and T^ 3 j(p 3 q 3 ) in equation O 
are coupled together. The number of coupled equations can be reduced since some 
of the particles are identical. In the case of three identical particles (NNN and AAA 



Table 12. Coupled channels (I, s) that contribute to a given NA or AA state 
with total angular momentum j and isospin i. 







NA 


AA 


3 


i 


M 


(l,s) 





1 




(0,0), (2,2) 


1 







(0,1), (2,1), (2,3) 


1 


1 


(0,1), (2,1), (2,2) 




1 


2 


(0,1), (2,1), (2,2) 


(0,1), (2,1), (2,3) 


2 


1 


(0,2),(2,1),(2,2) 


(0,2), (2,0), (2,2) 


2 


2 


(0,2), (2,1), (2,2) 







3 




(0,0), (2,2) 


3 







(0,3),(2,1),(2,3) 


2 


3 




(0,2), (2,0), (2,2) 


3 


2 




(0,3), (2,1), (2,3) 
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systems) all three amplitudes are equal and therefore equation l|7.3[l becomes, 

roc /> 1 

Ti-siiPili) = ^2 < S i i ^\ s 3 i 3 >si J Qjdqj J dcosBt^Xv^v'^E - qf /2^) 



Sjlj 

1 



T-mwi)- ( 7 - 4 ) 



Table El details the three NNN states characterized by total spin and isospin (S, I) 
that are possible as well as the two-body NN channels that contribute to each state. 
The 25 possible AAA states as well as the two-body AA channels that contribute to 
each state are also given. 

When two particles are identical and one different (NNA and NAA systems) 
two of the amplitudes are equal. With the assumption that particle 1 is the different 
one and particles 2 and 3 are the two identical, only the amplitudes ^^(piqx) and 
^•sffefe) are independent from each other. The reduction procedure for the case 
where one has two identical fermions has been described in [156] and [233] and will 
not be repeated here. Tables El an d El resume the possible NNA and NAA states 
characterized by its total spin and isospin (S, I) as well as the two-body NA and NN, 
or NA and AA, channels that contribute to each state. 

7.2.3. Results Before going to the numbers obtained two technical precisions are in 
order. It has been already discussed in section ITTI that, the (j, i) = (1,1) and (2,2) 



Table 13. Two-body AA (NN) channels (j, i) that contribute to a given AAA 
(NNN) state with total spin S and isospin /. 

NN AA 

s i CM) {hi) 

1/2 1/2 (1,0), (0,1) (1,2), (2,1) 

1/2 3/2 (0,1) (1,0), (1,2), (2,1), (2,3) 

3/2 1/2 (1,0) (0,1), (1,2), (2,1), (3,2) 

1/2 5/2 - (1,2), (2,1), (2,3) 

5/2 1/2 - (1,2), (2,1), (3,2) 

1/2 7/2 - (1,2), (2,3) 

7/2 1/2 - (2,1), (3,2) 

1/2 9/2 - (2,3) 

9/2 1/2 - (3,2) 

3/2 3/2 - (0,1), (0,3), (1,0), (1,2), (2,1), (2,3), (3,0), (3,2) 

3/2 5/2 - (0,1), (0,3), (1,2), (2,1), (2,3), (3,2) 

5/2 3/2 - (1,0), (1,2), (2,1), (2,3), (3,0), (3,2) 

3/2 7/2 - (0,3), (1,2), (2,3), (3,2) 

7/2 3/2 - (2,1), (2,3), (3,0), (3,2) 

3/2 9/2 - (0,3), (2,3) 

9/2 3/2 - (3,0), (3,2) 

5/2 5/2 - (1,2), (2,1), (2,3), (3,2) 

5/2 7/2 - (1,2), (2,3), (3,2) 

7/2 5/2 - (2,1), (2,3), (3,2) 

5/2 9/2 - (2,3) 

9/2 5/2 - (3,2) 

7/2 7/2 - (2,3), (3,2) 

7/2 9/2 - (2,3) 

9/2 7/2 - (3,2) 

9/2 9/2 - 
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Table 14. Two-body NA (jnAi^na) and NN channels (.Jnn>2nn) that 
contribute to a given NNA state with total spin S and isospin /. 



S I C?na,«na) (?nn,*nn) 



1/2 


1/2 


(1,1) 




1/2 


3/2 


(1,1),(1,2) 


(1,0) 


1/2 


5/2 


(1,2) 




3/2 


1/2 


(1,1),(2,1) 


(0,1) 


3/2 


3/2 


(1,1), (1,2), (2,1), (2,2) 


(1,0), (0,1) 


3/2 


5/2 


(1,2), (2,2) 


(0,1) 


5/2 


1/2 


(2,1) 




5/2 


3/2 


(2,1),(2,2) 


(1,0) 


5/2 


5/2 


(2,2) 





NA states, and the (j, i) = (2,3) and (3,2) AA states present quark Pauli blocking 
for L = 0, what translates into a strong repulsive core in the baryon-baryon potential. 
The reason for such a repulsive core is the presence of a forbidden state, that in the 
RGM potential should be eliminated from the relative motion wave function for each 
partial wave. This procedure is tedious both from the conceptual and numerical point 
of view [70,234]. It has been demonstrated [27] that for the Pauli blocked channels the 
local NA and AA potentials reproduce the qualitative behaviour of the RGM kernels 
after the subtraction of the forbidden states. This is why reference [225] used the local 
version of the quark Pauli blocked channels. 

For the three-body systems the binding-energy spectrum, say, the energy of the 
states measured with respect to the three-body threshold, as well as the separation- 
energy spectrum, the energy of the states measured with respect to the threshold of 
one-free particle and a bound state of the other two, will be calculated. The deepest 
bound three-body state is not the one with the largest binding energy but the one 
with the largest separation energy, since such a state is the one that requires more 
energy in order to become unbound, that is, to move it to the nearest threshold. The 



Table 15. Two-body NA (j'nAi'na) an d AA channels C?aa,*Aa) that 
contribute to a given NAA state with total spin S and isospin I. 



s 




O'na 


«na) 


0'aa, «aa) 


1/2 


1/2 


(1,1), 


(1,2), (2,1), (2,2) 


(1,0 


,(0,1) 


1/2 


3/2 


(1,1), 


(1,2), (2,1), (2,2) 


(0,1 


,(1,2) 


1/2 


5/2 


(1,1), 


(1,2), (2,1), (2,2) 


(0,3 


,(1,2) 


1/2 


7/2 


(1,2), 


(2,2) 


(0,3 




3/2 


1/2 


(1,1), 


(1,2), (2,1), (2,2) 


(1,0 


,(2,1) 


3/2 


3/2 


(1,1), 


(1,2), (2,1), (2,2) 


(1,2 


,(2,1) 


3/2 


5/2 


(1,1), 


(1,2), (2,1), (2,2) 


(1,2 


,(2,3) 


3/2 


7/2 


(1,2), 


(2,2) 


(2,3 




5/2 


1/2 


(1,1), 


(1,2), (2,1), (2, 2) 


(2,1 


,(3,0) 


5/2 


3/2 


(1,1), 


(1,2), (2,1), (2,2) 


(2,1 


,(3,2) 


5/2 


5/2 


(1,1), 


(1,2), (2,1), (2,2) 


(2,3 


,(3,2) 


5/2 


7/2 


(1,2), 


(2,2) 


(2,3 




7/2 


1/2 


(2,1), 


(2,2) 


(3,0 




7/2 


3/2 


(2,1), 


(2,2) 


(3,2 




7/2 


5/2 


(2,1), 


(2,2) 


(3,2 




7/2 


7/2 


(2,2) 
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numerical results for two-body and three-body systems appear in tables ED and IT7I 
respectively. 

7.2.4. The NN system Out of the two states of table [TT] only the (j, i) = (1, 0), that 
is the deuteron, is bound. The non-local model gives a deuteron binding energy of 
2.14 MeV, while the local version gives 3.13 MeV. The exact chiral constituent quark 
model NN potential gives a deuteron binding energy of 2.225 MeV [135]. This value 
was obtained by taking into account the A A partial wave (^aAj saa) = (4, 3) coupled 
together in addition to those given in table ^] The present calculation only considers 
S- and Z3-waves obtaining for the non-local potential a deuteron binding energy of 
2.14 MeV, which differs less than 0.1 MeV from the exact result. 

7.2.5. The NA system Out of the four possible NA states of table IT2l only one, 
the — (2, 1), has a bound state which lies exactly at the NA threshold for the 
local model, and with a small binding of 0.141 MeV for the non-local one. The states 
(ii *) = (1) 1) an d (2, 2) are unbound because they present quark Pauli blocking that, as 
in the case of the AA system, will play an important role in the three-body spectrum. 
The state (j, i) = (2, 1) can also exist in the NN system and there it corresponds to the 
1 I?2 partial wave which has a resonance at an invariant mass of 2.17 GeV [235-237]. 
This means that the NA bound state may decay into two nucleons and appear in the 
NN system as a resonance. As the NA bound state has for both local and non-local 
models energies very close to the NA threshold, the invariant mass of the system 
should be very close to 2.17 GeV. Thus, the chiral constituent quark model predicts 
the NN 1 £>2 resonance as being a NA bound state. 

7.2.6. The A A system Out of the eight possible A A states given in table ITSl with 
the non-local interaction five have a bound state, whereas the local interaction binds 
six of them (in both cases there are no excited states in any channel). The channels 
(j, i) = (2, 3) and (3, 2) are unbound as expected, due to the strong repulsive barrier 
at short distances in the S'-wave central interaction. These repulsive cores will largely 
determine the three-body spectrum. 

The predicted bound states = (1,0), (0,1), (2,1) and (3,0) appear also in the 
NN system. One can see that the deepest bound state is = (1,0), the second 



Table 16. Binding energies B 2 (in MeV) of the NN, NA and AA states, 'ub' 
denotes unbound channels. B^ (B^ L ) refers to the local (nonlocal) potential. 









NN 




NA 


AA 


j 


i 


n 2 


R NL 
D 2 


D 2 


R NL 
D 2 




a 2 





1 


ub 


ub 






108.4 


159.5 


1 





3.13 


2.14 






138.5 


190.3 


1 


1 






ub 


ub 






1 


2 






ub 


ub 


5.7 


ub 


2 


1 






0.0 


0.141 


30.5 


7.4 


2 


2 






ub 


ub 









3 










0.4 


0.2 


3 













29.9 


7.8 


2 


3 










ub 


ub 


3 


2 










ub 


ub 
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Figure 36. Real and imaginary parts of the the 3 Z?3 amplitude [236]. 

is (j, i) = (0,1), the third is (j, i) = (2,1), and the fourth is = (3,0). The first 
three states appear also, and precisely in the same order, in the NN system. The (j, i) 

— (1,0) state has the quantum numbers of the deuteron, the (j, i) — (0,1) is the 
virtual bound state, and the (j, i) — (2,1) is the 1 Z?2 resonance that lies at ~ 2.17 
GeV [235] (notice that the 3 i<3 NN resonance has no counterpart in table 1161 because 
only even parity states have been calculated and 3 F 3 has odd parity). Thus, the (j, i) 

— (3,0) state which is also allowed in the NN system could correspond to a new NN 
resonance found in reference [238]. 

It is interesting to note that some indication of a {j, i) = (3,0) resonance can 
already be seen in the recent analyses of the NN data by Arndt et al [236]. This 
channel corresponds to the NN 3 Ds partial wave. The most distinctive feature of a 
resonance is that as the energy increases the real part of the scattering amplitude 
changes sign from positive to negative while the imaginary part becomes large, so 
that the amplitude describes a counterclockwise loop in the Argand diagram. The 
energy at which this change of sign occurs corresponds to the mass of the resonance. 
Figure OBI shows the real and imaginary parts of the 3 Z)3 amplitude obtained from the 
single-energy analysis of reference [236]. As one can see a resonance-like behaviour 
seems to be present at about 700 and 1100 MeV. These kinetic energies correspond 
to invariant masses of 2.20 and 2.37 GeV so that in either case the ordering of the 
state agrees with that predicted in table El The A A bound state in the (j, i) = (3,0) 
channel has been predicted also by other models [79,95-98,239], and a method to 
search for it experimentally has been proposed [220]. 

7.2.7. The NNN system As a reliability test of the calculating scheme for the three- 
baryon system, the NNN bound-state problem has been studied. Of the three NNN 
states of table 1131 only the channel (S, I) = (1/2, 1/2), that is the triton, has a bound 
state. By using the local potential a binding energy of 5.76 MeV is obtained, and 6.52 
MeV with the non-local one. For comparison, we notice that the triton binding energy 
for the Reid-soft-core potential in the truncated T-matrix approximation is 6.58 MeV. 
Since the experimental value is B cxp = 8.49 MeV, the difference with the theoretical 
result, of about 3 MeV, is a measure of the uncertainty of the two-channel calculation 
in the case of the three-baryon system. A detailed study of the triton binding energy 
will be presented in the next section. 
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Table 17. Binding energies B 3 (in MeV) of the NNN, NAA, NNA, and AAA 
states, 'ub' denotes unbound channels. For the bound channels the separation 
energy is given below in parenthesis. (B^ L ) refers to the local (nonlocal) 

potential. 







NNN 




AT A A 

IN A A 




AT AT A 

IN IN A 


AAA 


s 


I 




tdN L 
B 3 








b? l 


B% 


B 3 


1/2 


1/2 


5.76 
(2.63) 


6.52 
(4.38) 


- 


- 


- 


- 


84.0 

(53.5) 


16.6 
(9.2) 


1/2 


3/2 


ub 


ub 


- 


- 


- 


- 


139.2 
(0.7) 


ub 


3/2 


1/2 


ub 


ub 


- 


- 


ub 


0.14 
(0.002) 


109.5 
(1.1) 


ub 


3/2 


3/2 






- 


- 


ub 


2.28 
(0.14) 






1/2 


5/2 






ub 


0.63 
(0.43) 










5/2 


1/2 






ub 


8.16 
(0.36) 






39.1 
(8.6) 


9.3 
(1.9) 


5/2 


5/2 






ub 


0.18 
(0.04) 










1/2 


7/2 














6.3 
(0.6) 


ub 


7/2 


1/2 














31.7 

(1.2) 


7.8 

(0.4) 


7/2 


3/2 














35.1 

(4.6) 


9.8 

(2.0) 



7.2.8. The NNA system Though the presence of bound states in the NA (j, i) — (2, 1) 
and NN (j, i) = (1,0) channels could induce to expect several bound states in the NAA 
system, this is not the case indeed. As a matter of fact, with the non-local potential 
only two out of nine possible three-body states given in table 1141 arc bound. Because 
of the attractive contribution of the NA (j, i) — (2, 1) bound state with the non-local 
model, the three-body state (3/2,1/2) turns out to be very weakly bound, with an 
energy of 0.14 MeV, and a separation energy different hardly from zero. This means 
that the (S, I) = (3/2,1/2) state is very near the NNA threshold and therefore it 
represents the tribaryon resonance with the lowest possible mass since it can decay 
into three nucleons and one pion. Also, for this case the three-body state (3/2,3/2) 
is bound. As it can be seen from table El this state has the contribution of all the 
two-body NA and NN channels. In spite of the fact that the NA two-body channels 

= (1, 1) and (2,2) present Pauli blocking, the attractive contribution of the NA 
Cm) = (2,1) and NN = (1,0) channels results to be enough to weakly bound 
this state. We note that neither one of the three-body states (S, I) — (3/2, 1/2) and 
(3/2,3/2) is bound with the local interaction. 

The predicted NNA states with (S,I) = (3/2, 1/2) and {S,I) = (3/2,3/2) which 
correspond to M k 3.4 GeV are the tribaryon resonances with the lowest mass and 
therefore the ones that would be more easy to detect experimentally. 
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7.2.9. The AAA system In the non-local case one finds that three of the sixteen 
possible NAA states given in tableED are bound. They are the (S, I) = (1/2,5/2), 
(5/2,1/2), and (5/2,5/2) states and their bound state energies are 0.63 MeV, 8.16 
MeV, and 0.18 MeV, respectively. In the case of the states (S,I) = (1/2,5/2) and 
(S, I) = (5/2,1/2) the repulsive barrier due to the quark Pauli blocked NA states 
(J, i) — (1, 1) and (2, 2) is less strong than the attraction due to the state (j, i) — (2, 1), 
so that they result to be bound in the non-local model. The state (S, /) = (5/2, 5/2) 
is the weakest bound state of this system, since in addition to the contribution of the 
NA quark Pauli blocked channels, there exists that of the AA quark Pauli blocked 
channels (J, i) — (2, 3) and (3, 2). This confirms what has been mentioned that it is the 
structure of the interaction of the two-body system the one which largely determines 
the three-body spectrum. Thus, the non-local interaction predicts the bound states 
(S,I) = (1/2,5/2), (5/2,1/2) and (5/2,5/2), which in principle may be observed as 
tribaryon resonances which decay into three nucleons and two pions with masses close 
to the NAA threshold. 

7.2.10. The AAA system While the non-local model presents four bound states the 
local interaction finds seven. From table 1171 one observes that the three states which 
are missing in the non-local version are barely bound in the local version, i.e., they 
have very small separation energies. Since the non-local interaction tends to lower the 
attraction in all the AAA channels it is not surprising that those which were barely 
bound have disappeared. The more strongly bound three-body state (that is, the one 
with the largest separation energy) is the (S, I) = (1/2, 1/2) state which has precisely 
the triton quantum numbers. This shows again, like in the AA and NN systems, the 
similarity between the AAA and NNN systems. 

The reason why the (S,I) = (1/2, 1/2) state is the more strongly bound is very 
simple. As shown in table^|this is the only state where none of the Pauli blocked two- 
body channels = (2,3) and (3,2) contribute. In all the other three-body states 
the strong repulsion of these channels either completely destroys the bound state or 
allows just a barely bound one. The state (S, I) = (7/2, 3/2) comes next with respect 
to the separation energy. This state has a somewhat anomalous behaviour since it 
has a relatively large separation energy. This behaviour is sort of accidental and it 
can be understood as follows. As seen in table there are four two-body channels 
contributing to the (S,I) = (7/2,3/2) state, the two attractive ones — (2,1) 

and (3,0) and the two repulsive ones = (2,3) and (3,2). However, as seen 

in table ITf)l the attractive channels (2,1) and (3,0) have bound states at E = —7.4 
MeV and E — —7.8 MeV, respectively, for the non-local model, and E — —30.5 
MeV and E = —29.9 MeV, respectively, for the local one, so that the poles in the 
scattering amplitude of these two channels are very close together and therefore there 
is a reinforcement between them, giving rise to the anomalously large separation energy 

7.3. The triton binding energy 

The relevance and/or necessity of considering the non-local part of NN potentials in 
realistic interactions is still under debate. Over the past few years several studies have 
appeared in the literature which stress the potential importance of non-local effects for 
the quantitative understanding of few-body observables and, specifically, for the triton 
binding energy [121,240-245]. However, the majority of these investigations [240-244] 
explore only non-localities arising from the baryonic meson-exchange picture of the 
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Table 18. Triton binding energy, Eg (i?Exp = —8.49 MeV), and wave function 
components: Pg is the L = symmetric spatial probability, Pgr the L = 
mixed symmetric spatial probability, Pp the L = 1 probability and Prj the L = 2 
probability. 





CCQM [80] 


Nijm II [248] 


Bonn B [249] 


£ s (MeV) 


-7.72 


-7.65 


-8.17 


Ps (%) 


91.49 


90.33 


91.35 


P S ' (%) 


1.430 


1.339 


1.368 


Pp (%) 


0.044 


0.064 


0.049 


Pd (%) 


7.033 


8.267 


7.235 



NN interaction. 

The non-localities generated in a quark-model derivation of baryonic potentials 
may be relevant for the three-nucleon bound state. It has been argued that the 
assumptions associated with meson-exchange models sharply limit the nature of 
the off-shell properties of those potentials, once the on-shell matrix elements are 
constrained to fit the NN data [246]. Therefore, it is very interesting to investigate 
the off-shell features of potentials derived from a quark-model. Some studies in this 
direction were done in reference [245] , but only the short-range part of the interaction 
was obtained by quark-model techniques, the intermediate and long-range parts being 
described by standard meson-exchanges between baryons. Accordingly, that work 
allowed only very limited conclusions with regard to effects of the quark substructure. 

The triton binding energy has been studied in detail by means of the chiral 
constituent quark model [80]. The calculation has been done in a momentum-space 
Faddeev formalism restricted to the 1 S'o and 3 Si — 3 Di NN partial waves, those which 
provide the bulk contribution to the three-nucleon binding energy. To solve the three- 
body Faddeev equations in momentum space a separable finite-rank expansion of 
the NN(coupled to NA) sector utilizing the EST method [247] is first performed. 
Such a technique has been extensively studied in other works for various realistic 
NN potentials [249] and specifically for a model including the coupling to the NA 
system [250]. In these works it was shown that, with a separable expansion of 
sufficiently high rank, reliable and accurate results on the three-body level can be 
achieved. In the present case it turns out that separable representations of rank 6 
(for 1 Sq — ( 5 Dq) and for 3 Si — 3 £>i) are sufficient to get converged results. We refer 
the reader to references [249-251] for technical details on the expansion method. The 
quality of the separable expansion on the NN sector is very good. The phase shifts 
for the original non-local potential and for the corresponding separable expansion are 
almost indistinguishable [80]. 

Results for the triton are summarized in table ^| It is reassuring to see that 
the predicted value of the triton binding energy is comparable to the values obtained 
from conventional NN potentials, such as the Bonn or Nijmegen models. Thus, these 
calculations show that quark model based NN interactions are definitely able to provide 
a realistic description of the triton. The results also give support to the use of such 
interaction model for further few-body calculations. One should not forget at this point 
that the number of free parameters is greatly reduced as compared, for example, to 
effective theories, and, in addition, they are strongly correlated by the requirement of 
describing the baryon spectrum. 

One can find in the literature other quark-inspired analysis of the triton binding 
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energy. Some preliminary studies were done in Ref. [245] with a hybrid model as 
mentioned above. Recently the triton has been studied with the same quark model 
approach of Ref. [47], used for the study of the NN phase shifts (see section l5.3.5[) . 
A 34-channel calculation is done. The result obtained, —8.390 MeV [121], is pretty 
close to experiment. Although this model is only intended for the study of the baryon- 
baryon interaction, i.e., the connection with the one-body problem is lost, the quality 
of its predictions is remarkable. 
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8. Summary 

We have reviewed recent calculations of few-body systems within a chiral constituent 
quark model framework. As for the moment an exact solution of QCD is not 
attainable, models incorporating its basic symmetries and some dynamics reveal as 
an important alternative to study low-energy hadron phenomenology In this energy 
regime, confinement and chiral symmetry breaking appear as the most important 
properties of QCD. When properly implemented by means of chiral constituent quark 
models, the interplay between gluons, Goldstone bosons and constituent quarks is 
able to provide a consistent description of the baryon structure and properties and 
the baryon-baryon phenomenology. Several aspects of chiral constituent quark models 
in their application to study few-body systems deserve to be emphasized. 

First, their universality in a double sense: on the one hand they allow to describe a 
great variety of phenomena (regarding N, NN and NNN systems) in terms of the same 
interacting potential, on the other hand once the very few parameters of the models 
are determined in the well-known NN system they allow to predict in a parameter-free 
way the spectrum and the interaction involving any other baryonic resonances. 

Second, quark antisymmetry is properly implemented showing up as cluster 
antisymmetry and quark exchange between clusters. Its interplay with quark dynamics 
provides an elegant explanation of the observed hard core of the NN force. Similarly a 
precise description of the quark Pauli short-range repulsion appearing in other systems 
is obtained. 

Third, these models at their current status have raised a degree of accuracy that 
allows to reproduce fine properties of two and three-baryon systems, like deuteron 
properties, the triton binding energy or charge symmetry breaking, with similar quality 
to baryonic meson-exchange models. Furthermore, the non-strange baryon spectrum 
is nicely reproduced once some rclativistic effects arc taken into account. 

Despite this success a lot of work is still to be done. The confining interaction 
lacks a fundamental description and the phenomenological constraints are not enough 
to determine its properties. This has an important influence on the spin-orbit force 
which for the moment is not correctly described. Moreover, a complete rclativistic 
treatment should be properly implemented and its consequences analyzed. Finally, 
other systems like mesons or multiquark systems should be studied within the same 
scheme. Meanwhile one can safely conclude that chiral constituent quark models 
have contributed in the last twenty years to a better understanding of the baryonic 
phenomenology. 
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